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Temperature and abundance profiles of hot gas in galaxy groups - 
II. Implications for feedback and ICM enrichment 
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ABSTRACT 

We investigate the history of galactic feedback and chemical enrichment within a sample of 
15 X-ray bright groups of galaxies, on the basis of the inferred Fe and Si distributions in the 
hot gas and the associated metal masses produced by core-collapse and type la supernovae 
(SN). Most of these cool-core groups show a central Fe and Si excess, which can be explained 
by prolonged enrichment by SN la and stellar winds in the central early-type galaxy alone, 
but with tentative evidence for additional processes contributing to core enrichment in hotter 
groups. Inferred metal mass-to-light ratios inside rsoo show a positive correlation with total 
group mass but are generally significantly lower than in clusters, due to a combination of lower 
global ICM abundances and gas-to-light ratios in groups. This metal deficiency is present 
for products from both SN la and SN II, and suggests that metals were either synthesized, 
released from galaxies, or retained within the ICM less efficiently in lower-mass systems. We 
explore possible causes, including variations in galaxy formation and metal release efficiency, 
cooling-out of metals, and gas and metal loss via AGN- or starburst-driven galactic winds 
from groups or their precursor filaments. Loss of enriched material from filaments coupled 
with post-collapse AGN feedback emerge as viable explanations, but we also find evidence 
for metals to have been released less efficiently from galaxies in cooler groups and for the ICM 
in these to appear chemically less evolved, possibly reflecting more extended star formation 
histories in less massive systems. Some implications for the hierarchical growth of clusters 
from groups are briefly discussed. 
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1 INTRODUCTION 

Groups of galaxies are repositories for a substantial fraction of 
all baryons and metals at low redshift. Being much more numer- 
ous than more massive clusters and more susceptible to the effects 
of non-gravitational processes such as galax y winds dri ven by su- 
pernovae or active galactic nuclei (e.g. . iHelsdon & PormiarjllOOOl : 
iBorgani et al.ll2004l ; lRomeo et al.ll2006h . they represent unique lab- 
oratories for the study of galactic feedback and chemical enrich- 
ment of the intergalactic medium. In practice, this can be exploited 
through X-ray observations of the hot intragroup gas, which allow 
the abundance of several heavy elements in this gas to be con- 
strained. 

Numerous such studies have been undertaken for mas- 
sive galaxy clusters, providing a reasonably well-constrained 
picture of the amount and distribution of metals in the in- 
traclus ter medium (ICM) and in the cluster galaxies them- 
selves jpinoguenov. David & Ponmanll200(]| : lB6hringer et al.l2004l : 
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De Grandi et all |2004 iMouhcine et al.l l2006l : ISaldi et alj l2007l : 
Leccardi & Molendil l20oir However, despite extensive observa- 



tional and theoretical efforts, several important questions related 
to ICM enrichment still remain, such as the relative importance of 
the mechanisms responsible for transferring metals from galaxies 
to the ICM, and the role of mergers and galactic wind activity in re- 
distributing enriched material. For example, it has been notoriously 
difficult to reproduce observed ICM metal distributions in cosmo- 
logical hydrodynamical simulations of clusters, suggesting that im- 
portant processes affecting baryons in these systems may still be in- 
adequatel y captured or altogether missing in such models (see, e.g., 
review bv lBorgani et al. 2008b). Characterization of the chemical 
makeup of the hot gas in more numerous lower-mass groups could 
provide an important ingredient in resolving these questions. Due 
to the structural role played by groups in connecting the 'field' - i.e. 
isolated galaxies and the intergalactic medium - with the rare high- 
density peaks represented by massive clusters, a clearer picture of 
enrichment in groups may have ramifications for our understanding 
of the evolution of baryons across a wide range of environments. 

In an early study o f the distribution of e l ement al abun- 
dances in intragroup gas, iFinoguenov & PonmanI ( Il999h investi- 
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gated three groups, for which they found significant spatial vari- 
ations in the ratio of iron to a-element abundances, conclud- 
ing that both supemovae (SNe) type la and II must have con- 
tributed to the enrichment of the intracluster medium in groups. 
Evidence for significant iron gradients in a sample of 10 groups 
and early-typ e galaxies were also found from RO SAT data by 
iBuotd ( l2000bh . lFinoguenov. Amaud & David ( 1200 ih followed this 
up using a larger sample of 18 cool clusters, but with only two 
systems at T < 2 keV, while iFinoguenov et alj feoOZ) studied 
nine groups mainly in the context of their entropy profiles and 
the constraints they imposed on preheating scenarios. More re- 
cent work include the two-di mensional XMM-Newton results of 
IFinoguenov et alj ( |2006| . |2007|) . but the main focus here was the en- 
tropy and pressure structure of the ICM. Thus, while most of these 
studies have focused on T < 2 keV systems, detailed discussions 
of the abundance distributions and the implications for galaxy and 
group evolution have generally been limited and based on samples 
of just a few groups (e.g.. IFinoguenov & Ponmanlll9 99). Hence, a 
dedicated study of chemical enrichment in a substantial number of 
groups has been lacking. 

In the first paper in this series jRasmussen & PonmanllZOOTl 
hereafter Paper I), we presented projected radial profiles of tem- 
perature and Fe and Si abundance of the hot intragroup gas in a 
sample of 15 groups, discussing the results in a statistical frame- 
w ork. The groups wer e selec ted mainly from the GEMS sample 
of 'Osmond & Ponman (2004), on the basis of their good photon 
statistics in available Chandra data and for their relatively undis- 
turbed X-ray morphology. The low Chandra background enabled 
robust abundance constraints to reasonably large radii, with the 
majority of our groups having Fe and Si profiles extracted to at 
least |r5oo, where rsoo is the radius enclosing a mean density of 
500 times the critical density. At the same time, the spatial reso- 
lution of Chandra allowed sampling of group cores at a level of 
detail usually unattainable with other X-ray telescopes. With only 
one exception (NGC4125), all the groups were found to contain 
cooler gas in the central regions. Main results of the study included 
the discovery that Fe abundances in the groups show a highly sig- 
nificant decrease with radius, declining continually from a central 
excess towards a value at rsoo which is, on average, lower than 
corresponding values seen in more massive clusters by a factor of 
two. In addition, it was found that Fe enrichment in group cores is 
dominated by SN la a s also found for more massive clusters (e.g., 
iDe Grandi etZll2004h. thus confirming early, tentative ASCA re- 
sults ( IFinoguenov et al.ll2000h . The SN II contribution to enrich- 
ment increases with radius, however, and potentially becomes ex- 
clusive at radii approaching rsoo ■ 

The purpose of the present paper is to explore some of the im- 
plications of these results. In Section |2| we outline the derivation 
of gas and metal masses and the optical properties of the groups. 
Section |3] presents the main results inferred from these quantities 
and from the observed abundance distributions, including the metal 
masses produced by SN la and SN II, the corresponding mass-to- 
light ratios, and the impact of SN in the brightest group galaxy 
(BGG) and in the groups as a whole. One finding is that the metal 
mass-to-light ratios deviate substantially from results for more mas- 
sive systems, and we explore a variety of possible explanations for 
this in Section|4| along with some further implications for the feed- 
back and chemical enrichment history of the groups in Section [5] 
Section |6| briefly discusses and summarizes the results, including 
some general implications for structure formation. 

As in Paper I, we use the solar abundance table of 
iGrevesse & Sauva theoretical supernova yields from 



llwamoto et"al] d 1 9991) and iNomoto et al] ( l2006ll . and the Kendall 
rank-order correlation coefficient r to quantify the significance ax 
of linear relationships. A value of Hq = 70 km s^^ Mpc^^ is 
assumed throughout, and all errors are given at the 68 per cent con- 
fidence level unless otherwise stated. 



2 METAL MASSES AND OPTICAL LUMINOSITIES 
WITHIN i?500 

Results from Paper I for the abundance profiles are summarised 
in Fig. [T] The data points here are similar to those of fig. 1 1 in 
Paper I, showing the results in radial bins of 20 measurements 
each. We briefly recapitulate the distinct features of these profiles, 
which are (i) a central abundance excess for both Fe and Si, of- 
ten accompanied by a dip in the very centre for both elements; 
(ii) a clear decline in Fe abundance with radius towards a value 
of ~ 0.1 Z0 at rsoo, with a somewhat less pronounced decline for 
Si (which furthermore flattens out or even rises at large radii). For 
reference, the SN yields adopted here, as explained in Paper I, are 
Mpe « 0.79 Mq and Msi « 0.21 M© for Fe and Si from SN la, 
and Mpc ^ 0.08 Mq and Afsi « 0.12 Mq for core-collapse super- 
novae including SN II. This implies Zsi/Z^c ~ 0.46 Zsi.e/ZpcQ 
for SN la ejecta, 2.62 Zsi.q/Zpc.q for those of SN II, and with a 
solar abundance ratio requiring 3.4 SN II per SN la. The impact of 
assuming different sets of SN yields will be discussed in SectionjS] 

2.1 Parametrizing the abundance profiles 

For the majority of our groups, abundances could not be measured 
all the way out to rsoo. In order to enable predictions of Z{r) 
beyond the radius of measurement for individual groups, as well 
as to facilitate visualisation of sample-averaged results, we made 
parametrizations of the average abundance profiles in Fig. [T] This 
approach can be justified on the grounds that while the statisti- 
cal uncertainties on Zyc and Zsi can be quite large for individual 
groups, the intrinsic scatter in these quantities is reasonably small 
outside the innermost regions. Outside the peak of the profiles in 
Fig.[T] we have adopted the fitted relations of Paper I [equations (8) 
and (9)] for these parametrizations. Inside the peaks, a linear rise 
in log Z vs. log (r/rsoo) has been assumed, normalised such as to 
ensure a continuous behaviour at the peaks. Outside r = O.lrsoo, 
where the vast majority of gas resides, the typical standard error 
on the mean in the individual bins in Fig.[T^ and[TJ) is 16 per cent 
for Zyc and 18 per cent for Zsi- We have conservatively adopted 
20 per cent relative errors as representative for our parametriza- 
tions of both profiles. These 20 per cent error envelopes are shown 
as shaded regions in Fig. [Tj within these errors, the parametriza- 
tions are consistent with every data point in the figure. Note that 
the shaded region in Fig. [T}; is not a direct parametrization of the 
Si/Fe ratios, but results from the ratio of the parametrized Si and Fe 
profiles, with their adopted errors added in quadrature. 

Unless otherwise stated, our analysis employs the observed 
abundance profiles of individual groups out to the edge of the us- 
able Chandra data, extrapolating these to larger radii where nec- 
essary using the above parametrizations and their associated er- 
rors. Though obviously a model-dependent approach, the adopted 
parametrizations are at least empirically motivated across all radii 
considered here, and the resulting Si/Fe ratios remain consistent 
with commonly used SN model yields everywhere inside rsoo ■ Fur- 
thermore, as will be demonstrated, there is no clear evidence for 
the observed abundance profiles at large radii to be systematically 
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Figure 1. Fe and Si abundance profiles for all groups, in radial bins of 20 data points, along with the corresponding Si/Fe ratios. Shaded areas illustrate the 
adopted Icr error envelopes on our parametrized abundance profiles. Dotted lines in (c) show the expectations for pure SN la and SN II enrichment. 



different for, for example, high- and low-temperature groups, so 
the adopted parametrizations should be reasonably representative 
for the entire group sample. The additional precaution of adopting 
conservative errors on these was taken to further alleviate concerns 
regarding this approach. 

2.2 Derivation of gas and metal masses 

For the calculation of gas and metal masses of the groups, we 
adopted a standard 1-D /3-model to describe the gas density 
profiles, using values of /3 and Tc for each group taken from 
the literature. These values are listed in Table Q] below, taken 
from the ROSAT analyses of Mulchaey et al. (2003) for NGC 507, 
NGC4125. and NGC7619 , lDavis et al.i ( ll996D for NGC 2300, and 
lOsmond & PonmanI ( [2004h for the remainder. 

The gas density profile of each group was normalised by 
means of the measured Chandra flux in the 0.1-0.3 rsoo range, 
using the spectral normalisation A from XSPEC, 

-|^g-14 /•0.3r5oo 

^= A.D-{l + .Y I ~dT/cm- (1) 

where D is the group distance (from Paper I) and and nu 
are the number densities of electrons and hydrogen atoms, respec- 
tively. The reason for normalising the profiles using the flux outside 
r = O.lrsoo is that many of the published gas density parameters 
result from two-component fits, with one component used to de- 
scribe the emission in and around the central galaxy, and the other, 
of interest here, to describe the extended group emission. This ap- 
proach also reduces the impact of projection effects on the result- 
ing normalisations, and suppresses effects related to blurring by 
the ROSAT po'ml spread function. Where relevant, a geometric cor- 
rection was applied to A, to account for lost detector area due to 
removed point sources and incomplete azimuthal coverage of the 
relevant radial range (see Paper I for details). The median correc- 
tion factor was 1.38, with no significant dependence on, for ex- 
ample, mean group temperature. Gas and metal masses were then 
calculated using the adopted /3-model parameters. Errors on these 
masses within the radius of interest were obtained on the basis of 
5,000 Monte Carlo realisations of the gas density and abundance 
profiles of each group. For each realisation, gas density parame- 
ters (A, P, Vc) and abundances (Zpc, Zsi) in each radial bin were 
drawn from Gaussians with a equal to the associated la error. The 
final metal mass and its uncertainty were then obtained as the mean 
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Figure 2. Illustration of the Monte Carlo output for the calculation of gas 
and metal masses and their uncertainties. Histograms show the obtained 
distribution of Fe masses for each group, with the resulting mean values 
represented by dashed lines, and the associated ilc uncertainty ranges by 
the shaded area. 

and (T of the resulting distribution of masses. These distributions 
are generally symmetric, as illustrated in Fig. [2] which shows an 
example outcome of this approach for the calculation of Fe masses. 

As is the case for the abundance profiles, the gas density pro- 
files have generally not been measured out to rsoo within our sam- 
ple. Where relevant, the adopted /3-parameters were therefore used 
to extrapolate the gas density profile of each group beyond the ra- 
dius of measurement. Systematic uncertainties related to this ex- 
trapolation are not straightforward to quantify, because the num- 
ber of gas mass measurements extending to rsoo in T < 2 keV 
systems are still fairly limited (see, e.g., ISun et al1l2009h. XMM- 
Newt on studies of groups to large radii feasmussen & PonmanI 
l2004h suggest that standard /3-models can provide a good descrip- 
tion of the profiles of relaxed groups out to, and beyond, rsoo . How- 
ever, several studies have found the slope of the gas density pro- 
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file in more massive clusters to steepen progressively at large radii, 
with typical slopes at rsoo that are '-^ 15 pe r cent steeper than the 
'can onical' cluster value of /3 ~ 2/3 (e.g.. Ivikhhnin et al ] |2006l : 
iMaug han et al. 2008). To account for this steepening, Vikhlinin 
et al. (2006; see also ISunetal.|[2()0^ multiply the standard /3- 
model for the density profile by a factor [(1 + r'' /r'JY^'']~^ , with 
7 = 3 and and e fitted parameters. 

Although the Vikhlinin et al. ( 2006) study was largely limited 
to (T) > 2 keV systems (leaving it unclear to what extent these re- 
sults can be reliably extended to the regime of cool groups, whose 
density profiles inside rsoo tend to be flatter than those of clusters; 
[Sanderson et al. 2003), we have investigated the impact of adopting 
the ^Vikhlinin et al. ( 2006) parametrization for our groups. Mean 
values for the Vikhlinin clusters are Ts — l.lrsoo and e = 3.2 
(from their table 2), with no systematic dependence of these pa- 
rameters on cluster temperature; adopting these values, we find that 
our gas masses within rsoo could, on average, be overestimated by 
24 per cent (with a small standard deviation of 4 per cent). This 
is within the adopted systematic uncertainty on A/gas (see Sec- 
tion |2!4ll, so within these adopted errors, our gas masses should 
be robust towards such changes in the assumed shape of the den- 
sity profile. Furthermore, any trends in A/gas-based quantities with 
mean system temperature would not be significantly affected by 
such a change. Many of the cluster comparison results presented 
in the following may also be subject to similar corrections, so 
our overall conclusions should remain valid irrespective of such 
a change. 



2.3 Optical, near-infrared, and radio properties 

In order to derive mass-to-light ratios and to estimate the expected 
supernova rate in the groups, K- and B-band luminosities of the 
group members were also derived inside rsoo. Luminosities were 
computed from extinction-corrected galaxy magnitudes in NED us- 
ing the method of Osmond & Ponman (2004): Galaxies were ex- 
tracted inside rsoo and within a velocity interval of 3(Tv around 
the mean recession velocity Vr of the group, using starting values 

Jsmond i 



of Vr and velocity disper sion g-y from 



Osmond & Ponmai 



ig values 
ri i2004l) 



where available, or from Mulchaey et al. (2003) otherwise. These 
quantities were then re-evaluated iteratively until convergence was 
reached. The only exception to this rule was NGC 1407, for which 
a velocity range extending to S.Gctv was employed. This was done 
in order to include the second brightest group galaxy, NGC 1400, 
which has an unusually lar ge peculiar velo city of ~ 1100 km s^^ 
relative to the group mean dTrentham. Tullv & Mahdavi 2006 ). 

The adopted isT-band magnitudes of the group members gen- 
erally derive f rom Two-Micron All Sky Survey data (2MASS; 
ISkrutskie et al. l 2006). Inspection of a composite histogram of the 
apparent magnitudes suggests that incompleteness sets in around 
rriK ~ 12.5 for our groups, consistent with the expectation that 
the 2MASS source catalog is 90 per cent complete down to rriK = 
13.5 for extended sources more than 30° from the Galactic plane 
dSkr utskie et al. 2006), a region covering 12 of our 15 groups. For 
each group, we corrected for incompleteness below the absolute 
magnitude correspond i ng to rriK = 12.5 using the Schechter fit of 
iLin. Mohr & Stanford ( |2004|) to the 7^-band luminosity function of 
93 groups and clusters (and including their correction factor to ob- 
tain the Schechter normalisation inside rsoo instead of r2oo)- The 
resulting corrections are typically small, on average 5 per cent. 

In order to facilitate comparison to existing results in the lit- 
erature, B-band magnitudes of the group members were also ex- 
tracted from NED in a similar manner. In four of the groups, mag- 



nitudes had to be converted from the _R-band for a few galaxies. 
This was done assuming B-R = 1.5 as is appropriate for early- 
type galaxies. While this could lead to a slight underestimation of 
LB,v/e found that using instead B — _R = 0.8 as for spirals resulted 
in a mere 4 per cent increase in Lb for NGC 4325, with the change 
being less than 0.5 per cent for the remaining three groups. We have 
refrained from taking these minute corrections into account in the 
error budget. Since the optical magnitudes generally derive from a 
variety of heterogeneous data sets with varying completeness, we 
have not attempted a similar correcti on for incompletenes s as ap - 
plied to the if -band data. Following lOsmond & Ponmai] ( l2004h . 
we instead simply aimed to homogenize completeness across the 
sample by applying an absolute magnitude cut to all groups. This 
cut discarded galaxies fainter than Mb = —16.32, retaining an 
expected 90 per cent of the B-band luminosity of the galaxy popu- 
lation. 

Finally, using the derived B- and A'-band luminosities and the 
corresponding B-K colours, estimates of the total stellar masses 
A/* of the groups and their central galaxies were obtaine d using 
the coefficients listed in table 1 of iBell & de Jond d200lh . which 
imply 



log 



M,/Lk 
M0/L0 



= 0.212(5 ~K)- 0.959 



(2) 



for the A'-band stellar mass-to-light ratio of the group members. 
We note that our approach implies that our B-K colours could be 
slightly but systematically biased towards the red, because we are 
attempting to account for all Jf-band light by performing a com- 
pleteness correction while only retaining an expected 90 per cent of 
the B-band light. This could lead to overestimates of stellar mass 
which, while potentially varying from group to group, should not 
exceed 0.02 dex, i.e. 5 per cent, and so this bias has no substantive 
impact on our results. 

2.4 Robustness of results 

Although not themselves the focus of the present study, reliable gas 
masses are still important ingredients in this work, so we briefly dis- 
cuss the robustness of our adopted approach. For example, while 
the derived statistical errors from our Monte Carlo calculations 
are small for most groups, one concern is whether the adopted /3- 
parameters and our method of normalising the density profile are 
robust. Another is that our approach does not take into account ra- 
dial variations of the plasma emissivity. To explore the potential im- 
pact of this, we have compared our gas mass estimates to published 
values for the groups for which we have been able to find useful 
comparison data in the literature, consisting of either Chandra, 
XMM, or Suzaku measurem ents published within a given radius 
using different approaches (iMorita et al.l |2006| ; iGastaldello et al.l 
|2007| ; iKomivama et alj|2009h . At larger radii, we have also com- 
pared our results within rspp to those ob t ained using the density 
parametrizations of lSanderson & PonmanI ( |2003|) based on ROSAT 
data. This comparison is shown in Fig.|3] 

The figure indicates that there is no clear systematic vari- 
ation between our results and other estimates, with most esti- 
mates agreeing to within ^ 25 per cent. An exception is the 
clearly discrepant NGC 5846, for which the lSanderson & PonmanI 
(2003) parametrization yields a gas mass within rsoo which is only 
65 per cent of our estimate. However, even when including this sys- 
tem, the average deviation of our estimates from existing results is 
only 10 per cent for all measurements, or 17 per cent if only in- 
cluding the independent largest-radius estimate for each group. We 
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Figure 3. Ratio of our estimated gas masses to other published results. 
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Figure 4. Comparison of (top) Fe and (bottom) Si profiles derived for 
NGC 4325 using the most recent Chandra calibration data (black crosses) 
and those of our original analysis in Paper I (gray diamonds). 



therefore consider a 25 per cent error a reasonably conservative es- 
timate of the systematic uncertainty associated with our method of 
gas mass estimation, so we have added this uncertainty in quadra- 
ture to those related to any quantity based on A/gas • 

Since the time of our original Chandra analysis, performed 
using CIAO 3.0 with calibration database CALDB 3.2.3, updated 
calibration data including a recalibration of the Chandra effective 
area have been released. To briefly test the potential impact of this 
update on our results, we rederived the temperature and abundance 
profiles for NGC 4325, a group typical of our sample in terms of 
signal-to-noise ratio and mean ICM temperature. For this updated 
analysis, we employed the most recent calibration data, CALDB 
4.1.2, with CIAO 4.1, following the spectral fitting procedure de- 
scribed in Paper I. We find that the updated temperature measure- 
ments are perfectly consistent with the original values in all radial 
bins, and the results for the Fe and Si profiles, shown in Fig. [4] 
also demonstrate good consistency throughout, with no systematic 
offset between the new and original results. In fact, the main effect 
of using the updated calibration data is a considerable reduction 
in the statistical uncertainties of most measurements. We therefore 
conclude that our original results should generally be robust to the 
recent changes to the Chandra calibration. 

As a further indicator of the robustness of our abundance mea- 
surements, we point to the good overall agreement wit h the in- 
dependently obtained XMM results presented by Finogue nov et al.l 
hood boOTi) (see discussion in Paper I), in particular the excel- 
lent agreement across the full radial range between our sample- 
averaged Fe profile in Fig. [T] and the corresponding XMM result 
(Johnson et al., in preparation). Encouragingly, good agreement 
is also seen when comparing our metal masses derived within 
specific radii to the corresponding XMM results for NGC 5044 
teuote, Brighenti"& Mathews! [20041) and the Suzaku results for 
HCG 62 i F^koi et alj l2008l). NG C 507 jSato et alj l2009t) . and 
NGC 5044 jKomivama et al.l2009l) . once differences in the adopted 
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Figure 5. Comparison of (left) our Fe masses and (right) mass-to-light ra- 
tios to existing results, where available. Dotted lines represent equality. 



abundance table have been taken into account. This is illustrated in 
Fig.[5](left panel), demonstrating consistency at the ~ la level be- 
tween ours and existing results in all four cases. Comparing also our 
resulting Fe mass-to-light ratios to the corresponding Suzaku esti- 
mates in Fig.[5](right panel), good overall agreement is again seen, 
with HCG 62 as the only exception. Since the derived Fe masses for 
this group are consistent, this discrepancy stems from differences 
in the assumed value of Lb within r = 13 arcmin, with the value 
of lTokoi etal]j2008l) . Lb = 4.9 x 10^" (derived by converting 
i?-band magnitudes to the B-band assuming a fixed B-R colour), 
being significantly lower than our value of Lb ~ 7.6 x 10^" 
obtained with the approach outlined in the previous Section. This 
minor issue notwithstanding, we conclude that our results for both 
gas masses, metal masses, and metal mass-to-light ratios generally 
show good agreement with existing results in the literature, thus 
lending credibility to our approach. 
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Table 1. Gas density parameters, gas and metal masses, optical and near-infrared luminosities, and galaxy velocity dispersions for the sample, all obtained 
inside rsoo- Also listed is the corresponding fractional metal mass fi^ provided by SN la for both Fe and Si. 



Group 




P 






-^-^gas 




/la 




/la 


Lb 


Lk 


0"v 




(keV) 




(kpc) 




(IOI2 Mq) 


(10« Mq) 


(Fe) 


(10* Mq) 


(Si) 


(logLo) 


(log L0) 


(km s-l) 


NGC383 


1.65 A nt^ 
— U.UD 


0.36 ± 0.00 


2.1 ±0, 


.2 


5.4 ±0.1 


18.4 ±4.7 


0.35 


21.2 ± 7.9 


0.08 


11.50 


12.26 


450 


NGC507 


1-30+^3 

— U.Uo 


0.42 ± 0.01 


25.7 ±0, 


,6 


3.9 ± 0.2 


14.4 ±3.9 


0.31 


16.4 ±6.3 


0.07 


11.69 


12.35 


635 


NGC533 


1 99+0.05 


0.42 ± 0.01 


2.2 ± 1, 


.7 


2.3 ± 0.1 


5.0 ± 1.4 


0.25 


6.2 ± 2.7 


0.06 


11.47 


12.11 


439 


NGC741 


1 42+0.14 


0.44 ±0.01 


2.3 ±0, 


.2 


2.7 ±0.1 


3.5 ± 1.3 


0.34 


11.9 ± 3.8 


0.03 


11.35 


12.03 


453 


NGC 1407 


1 01+"°'' 


0.46 ±0.01 


0.1 ±0, 


.1 


0.7± 0.1 


1.7±0.3 


0.31 


2.0 ±0.6 


0.07 


11.04 


11.73 


319 


NGC2300 




0.41 ±0.03 


56.5 ± 16, 


.8 


1.2 ± 0.4 


2.8 ± 1.0 


0.19 


3.8 ± 1.8 


0.04 


10.86 


11.49 


300 


HCG42 


o.solH^ 


0.56 ±0.02 


4.7 ±0, 


.7 


0.6 ±0.1 


1.2 ±0.4 


0.38 


1.3 ± 0.6 


0.10 


11.32 


11.92 


282 


MKW4 


1 70+0.07 


0.43 ±0.01 


9.4 ± 2, 


.9 


7.0 ±0.2 


22.3 ±6.6 


0.54 


19.7 ±9.7 


0.16 


11.78 


12.44 


565 


NGC 4125 




0.47 ±0.03 


< 0, 


.7 


0.2 ± 0.1 


0.4 ±0.1 


0.27 


0.6 ±0.3 


0.05 


10.62 


11.29 


56 


NGC 4325 


qg+0.02 


0.58 ±0.01 


27.6 ±5, 


.0 


1.8 ±0.1 


4.3 ± 1.4 


0.71 


2.8 ± 1.5 


0.30 


10.96 


11.73 


376 


HCG62 


^•'-"-'-0.03 


0.48 ±0.01 


2.4 ±0, 


.3 


3.1 ± 0.1 


5.3 ± 1.4 


0.52 


10.7 ±3.5 


0.06 


11.35 


12.02 


418 


NGC 5044 


1 19+0.03 
1-^^-0.03 


0.51 ±0.00 


6.0 ±0, 


.2 


2.3 ±0.1 


6.3 ± 1.5 


0.42 


6.6 ±2.0 


0.10 


11.22 


11.92 


426 


NGC 5846 


ee+'^o" 

U-DD_0.03 


0.51 ±0.01 


2.2 ±0, 


.3 


0.9 ±0.1 


2.1 ±0.4 


0.29 


2.6 ±0.8 


0.06 


11.23 


11.90 


346 


NGC 6338 


9 1 n+0.19 

^•-^'-'.-0.07 


0.44 ±0.01 


3.7± 1, 


.0 


9.4 ±0.3 


22.6 ±10.2 


0.55 


34.6 ± 14.7 


0.05 


11.54 


12.37 


651 


NGC 7619 


^•^"-0.03 


0.44 ±0.01 


7.9 ±0, 


.5 


1.5 ± 0.1 


3.2 ±0.7 


0.45 


3.2 ± 1.2 


0.12 


11.44 


12.12 


557 



3 RESULTS 

This section summarizes the derived group properties that will form 
the basis of the discussion in the following sections. Table[T]lists the 
main group properties used in this study, including the mean tem- 
perature (T) of each group derived inside the range 0.1 — O.Srsoo 
(from Paper I), the adopted gas density parameters, derived metal 
masses inside rsoo and their fractional contribution from SN la, the 
total B- and A"-band luminosities inside the same radius, and the 
derived galaxy velocity dispersions. Some of the properties of the 
brightest group galaxy (BGG) in each system (in all cases an early- 
type galaxy at the centre of the X-ray emission) are also relevant to 
the following discussion. In addition to stellar masses derived using 
equation we have also compiled central stellar velocity disper- 
sions ao,. for these galaxies from Hyperleda (available for all the 
BGGs except NGC 4325 and 6338), and the 1.4-GHz luminosity 
density of any central radio source, based on fluxes listed in NED. 
These results are listed in Table |2] 

Based on the parametrized Si/Fe ratios in Fig. [T] and the 
adopted SN yields, the abundance profiles have been decomposed 
into sample-averaged contributions from the two SN types in Fig.|6] 
Errors in this diagram were propagated from the adopted statistical 
errors on our parametrizations of Z-pcif) and Za,i(r) but are domi- 
nated at most radii by those related to the SN yields. The latter were 
evaluated by varying the assumed SN yields accord ing to the dif - 
fere nt SN la and SN II models tabulated by Iwamot o et alj h999l) 
and iGibson. Loewenstein & Mushotzkvl l |l997,) , respectively, sub- 
ject to the constraint Zsi/Z-pe > 2.0 for SN II, as required by our 
results at large radii in Fig. [T] At each radius, the resulting max- 
imum deviation relative to the case of our "reference" yields was 
taken as the uncertainty associated with the assumed SN yields and 
then added in quadrature to the statistical errors on the abundance 
parametrizations. The shaded regions in Fig.[6]outline the resulting 
lower and upper limits. The large uncertainties for Si in particular 
derive mainly from the significant variations in predicted SN la Si 
yields (0.14-0.27 Mq) among the models considered. Nevertheless, 
at most radii r > O.lrsoo considered in this study, the resulting un- 
certainties can be accounted for by assuming a characteristic error 
of 35 per cent, illustrated by the fiducial error bars in Fig.|6] Wher- 



Table 2. Summary of relevant properties of the brightest group galaxy, in- 
cluding the derived stellar mass, central stellar velocity dispersion, and 1 .4- 
GHz radio power of any central radio source. 



BGG Lk M* o-q,, P1.4 

(logLo) (log Mo) (kms-i) (WHz-i) 



NGC 383 


11.66 


11 


.69 


278 


3.4 


X 


10- 


NGC 507 


11.69 


11 


.55 


308 


3.7 


X 


10- 


NGC 533 


11.72 


11 


.56 


273 


2.2 


X 


10- 


NGC 741 


11.81 


11 


.65 


295 


7.2 


X 


10- 


NGC 1407 


11.48 


11 


.37 


273 


6.8 


X 


10- 


NGC 2300 


11.24 


11 


.15 


261 


4.8 


X 


10- 


HCG42 


11.71 


11 


.57 


320 








MKW4 


11.90 


11 


.77 


275 


1.9 


X 


10- 


NGC 4 125 


11.27 


11 


.12 


227 


1.7 


X 


10- 


NGC 4325 


11.35 


11 


.26 










HCG62 


11.62 


11 


.65 


232 


1.6 


X 


10- 


NGC 5044 


11.33 


11 


.13 


237 


5.0 


X 


10- 


NGC 5846 


11.51 


11 


.39 


237 


2.3 


X 


10- 


NGC 6338 


11.78 


11 


.67 




1.1 


X 


10- 


NGC 7619 


11.55 


11 


.42 


322 


7.3 


X 


10- 



ever relevant, we have taken this uncertainty into account in the er- 
ror budget by adding it in quadrature to existing errors, for example 
when computing total metal masses contributed by either SN type. 
We also stress that while Fig.[6]remains useful as a representation of 
the sample-averaged results, the SN decomposition for individual 
groups was performed using the actual measurements pertaining 
to each group, with the results in the Figure being employed only 
when extrapolation to 7-500 was necessary. 

In accordance with res ults for clusters, Fi g. [6| establishes on 
a firmer basis the results of iFinoguenov et alj [2000 .) that SN la 
products in groups are centrally concentrated, contributing about 
80 per cent of the Fe in group cores, and that SN II products on 
average are more uniformly distributed, showing only mildly in- 
creased central levels. Fe produced by SN II therefore dominates at 
large radii both in groups and clusters. Fe enrichment levels in the 
central metal excesses are comparable to those in cool-core clus- 
ters and are clearly dominated by SN la as also found in clusters 
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Figure 6. Sample-averaged contributions to the total abundance of (a) Fe 
and (b) Si from SN la and SN II, based on our Z{t) parametrizations. EiTor 
bars at O.lrsoo and at large radii show the adopted characteristic uncertain- 
ties of 35 per cent. 

jPe Grandi et alj|2004lde Plaa et alj|2007l) . but cannot be entirely 
explained by a central excess in SN la products alone unless assum- 
in g non-standard SN la yields. T his agrees with the earlier result 
of iFinoguenov & PonmanI ( 1 19991) based on a sample of only three 
groups, as does the derived SN II contribution to the Fe abundance 
at large radii of 0.05-0. 1 Zq . 

Total gas masses normalised by optical luminosity are sum- 
marized in Fig. [7] The figure reveals large variations in the ratio of 
gas mass to B-band light, and suggests a systematic dependence on 
(T), albeit one subject to significant scatter (Kendall's r = 0.35, 
correlation significance ok ~ +1.8, which reduces to ax = 1-6 if 
using Lk instead of Lb)- Inferred values range from 3~30 Mq/Lq, 
with a mean of 12 ± 2 Mq/Lq for groups in the T = 1 ~ 2.2 keV 
range. Although the latter val ue is in good agreem ent with results 
of cosmological simulations jRomeo et alj|200q) . it is nominally 
somewh at lower than the corre sponding v alue of ~ 20 Mq/Lq 
found bv IFinoguenov. Burkert & Bohringeil ( l2003h . and the overall 
implication of Fig.|7]is that gas mass-to-light ratios in our groups 
are well below typical cluster values. 



3.1 Metal masses and mass-to-light ratios 

The results for the total X-ray metal masses in the groups are sum- 
marized in Fig. [8] including the decomposition into contributions 
from SN la and SN II. Also shown are the corresponding metal 
mass-to-light ratios, all presented as a function of mean group tem- 
perature (T). We note that (T) correlates with total Lb, Lk, M,, 
and velocity dispersion at 3.6(7-3. 8a significance, suggesting 
that none of these parameters is a seriously biased proxy for total 
group mass. Also note that we generally normalise by Lb rather 
than Lk or stellar mass in these plots in order to enable straight- 
forward comparison to existing results in the literature. Since total 
Lk and Lb are tightly correlated for the sample (ok ~ 4.9), this 




<T) (keV) 

Figure 7. Gas mass-to-light ratio i nside rspp . Dashed line s hows the mean 
value of the T > 4 keV clusters of lFinoguenov et ei\ )2003l) . 



particular choice should have little bearing on the qualitative trends 
or the degree of intrinsic scatter seen in the figures. 

As can be seen from Table [T] and Fig. [8] the total Fe and Si 
masses derived inside rsoo span two orders of magnitude across 
the sample, roug hly from M ~ 3 x 10''-3 x 10^ Mq for both el- 
ements, and with a clear dependence on total group 'mass' (T) as 
anticipated. For Si, the higher abundance observed at large radii (cf. 
Fig. [T} roughly compensates for its lower atomic mass relative to 
Fe, yielding rather similar aggregate masses for the two elements 
in most cases. The fractional contribution by SN la to the metal 
production as listed in Table[T]implies that, on average, SN II have 
generated ^ 60 per cent of the Fe and ~ 90 per cent of the Si in- 
side rsoo. There is a weak tendency for these fractions to decrease 
with (T), (jk = —1.6 (aK ~ —2.1 if excluding the clearly dis- 
crepant NGC 4325), suggesting that enrichment in cooler systems 
is slightly more SN Il-dominated. 

The iron mass-to-light ratio (IMLR) represents a key quantity 
for studies of enrichment in groups and clusters, as it is expected to 
reflect the ability of a system to generate and retain enriched ma- 
terial within its gravitational potential. As already mentioned, the 
derived Fe abundances of our groups at large radii, where most of 
the gas resides, are lower than the corresponding cluster values by 
a factor of ~ 2. As discussed in Paper I, a similar deficit, if slightly 
smaller, is also observed for Si. The combination of a lower metal 
content in the hot gas and generally lower gas mass-to-light ratios 
in groups comp ared to more massive clusters (cf. Fig. |7] see also 
iLin et '^ 120031) should manifest itself in lower total metal M/L 
ratios in cool systems. For iron, we note th at a similar conclusion 
was already put forward bv iRenzinil ( 1 1 997h on the basis of ROSAT 
and ASCA data collected from the literature. However, this study 
used global abundance measurements out to the radius of X-ray 
detection, thus neglecting the important fact that X-ray emission in 
groups is generally detected to relatively smaller radii than clusters. 
Since cumulative Afpo/Ls ratios tend to rise outwards because the 
gas g enerally has a flatter distribution than tha t of the optical light 
(e.g.. IFinoguenov et al.ll2000l : ISato et al.ll2009h . the bias thus intro- 
duced will result in lower values in X-ray faint systems, even if 
the intrinsic Fe KI / L ratio inside a fixed fraction of the virial ra- 
dius did not vary between groups and clusters. The problem is only 
aggravated by the fact that groups typically have shallower gas dis- 
tributions than clusters and so contain a relatively larger fraction of 
their hot gas (and metals) at large radii. A fair comparison between 
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Figure 8. Left panel: Derived total masses of Fe and Si inside rsoo as a function of mean group temperature (T), along with the corresponding _B-band 
mass-to-light ratios. Centre and right panels show the results decomposed into contributions from SN la and SN II, respectively. Dashed lines show the mean 
values for the T > 4 keV clusters of .Finoguenov et al.. l .2Q03i) . and shaded regions mark the typical cluster ranges from .Finoguenov et aL ( .2OQ0i) . 



systems of differing X-ray brightness therefore requires all systems 
to be consistently compared within a fixed overdensity radius. 

These issues prompted us to investigate aggregate metal mass- 
to-light ratios of our groups for comparison to clusters. The results 
within rsoo are also shown in Fig. [8] The inferred ratios vary by a 
factor of ~ 10 across the sample, with indications of both AIpa/L 
(aK = 2.1) and Msi/L (ax = 2.3) systematically increasing 
with group temperature. Note that this is despite the fact that our 
temperature-independent Z{r) parametrizations, adopted for the 
groups for which the abundance measurements do not extend to 
rsoo, will tend to suppress trends with (T) (the main motivation 
for showing Fig [8] in its current form is to facilitate comparison 
to more massive systems, for which results are typically provided 
within rsoo)- The results clearly suggest lower Fe (and Si) mass- 
to-light ratios in cooler systems, with both ratios only approach- 



ing typical cluster values at the high-T end covered by our sam- 
ple. We note that this apparent shortfall of metals in groups is not 
simply an artefact of our choice of cluster comparison data . For 
example, the T > 3 keV clusters of lOe Grandi et all ( |2004) dis- 
play an IMLR within rsoo of 1.7 ± 0.3 x 10"^ for our adopted 
value of Hp, even high er than the ~ 1.3 x 10~^ reported by 
iFinog uenov et al.'('2003'). Using a considerably larger group sample 
than that of Finoguenov et al. (2000), we thus confirm their tenta- 
tive result that groups generally have Fe M/L of < 0.006 at all 
radii r < rsoo and tend to show lower Mpc/is than clusters at a 
fixed overdensity radius. 

In order to infer the origin of the lower metal mass-to-light 
ratios in cool systems, a key question to address is whether the 
shortfall of Fe compared to clusters is present both for SN la 
and SN II products. In Fig. [8]i and i, the IMLR has been decom- 
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Figure 9. Stacked Fe and Si profiles for 'cool' ((T) 5C 1 . 1 keV) and 'warm' 
((T) > 1.1 keV) subsamples separately, in radial bins of 10 data points. 



posed into contributions from SN la and SN II (denoted IMLRia 
and IMLRii in the following). This clearly reveals a deficit of 
Fe from both SN types in gr oups, a conclusion also reached by 
iFinoguenov & Porrniarj jl999l) , possibly s uggesting a combination 
of re sponsible mechanisms. We no te that IFinoguenov & Poimiarj 
( [1999,) and Finoguenov et'al] ( l200lh see a clear increase in IMLRn 
with T for clusters (T > 2 keV), and a global deficit of SN II ejecta 
in groups (as do we), but their cluster IMLRia is nearly constant 
with T across the full range covered by their systems, at variance 
with our results in the group regime. 



3.2 Dependence on ICM temperature 

Figs.|7]and[8]suggest that the observed variation in IMLR among 
the groups is driven primarily by a systematically increasing gas 
mass-to-light ratio with mean group temperature. In order to also 
investigate any temperature-dependence on abundance structure, 
the binned abundance profiles of Fig. [T] were split into results for 
'cool' (T s; 1.1 keV) and 'warm' (T > 1.1 keV) groups in Fig.[9] 
We repeat here that (T) was derived within the fixed radial range 
0.1 — O.Srsoo and so is not affected by the strong central cooling 
seen in most of the groups (cf. Paper I). Evidently, Fe abundances in 
cooler groups tend to be lower at all radii common for the two sub- 
samples. This conclusion is clearly robust inside r < O.Srsoo, as 
will be demonstrated quantitatively below, but the situation is less 
clear at larger radii where only a few of the 'cool' groups within 
our sample have useful Chandra data. The results for Si are qualita- 
tively similar to those of Fe in the group cores (r < 0.1 — 0.2r5oo), 
but at larger radii there are no obvious differences between the two 
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Figure 10. Mass-weighted Fe abundances within various radial ranges. Our 
results are marked with filled circles, with the results of IFinoguenov et aU 
12006, 2007) included as empty circles where relevant. Only groups with 
data covering the full radial interval are shown in each plot. Labels in lower 
right corners specify the significance afc of any linear correlation among 
our own data points. 



subsamples. Note also that the poorer statistics resulting from di- 
viding the sample masks the clear indication of a central abundance 
dip seen in Fig.[T] 

Fig. |9| indicates a general trend of higher metal abundances 
in hotter groups, at least in their central regions. To explore this 
quantitatively while also accounting for any systematic differences 
in the amount and distribution of gas with system temperature, we 
consider in Fig.[TO|the Fe mass-to-gas ratios, effectively the mass- 
weighted Fe abundances. The top panel shows the ratios within 
0.5^500, for which we have measurements for all groups and do not 
have to rely upon the (temperature-independent) radial abundance 
parametrizations. Within this radius, there is evidence for a system- 
atic increase in mean Zfc with (T) at a moderate significance level 
of 2.4-a. The following panels break this down into separate radial 
intervals, only including groups with data in the relevant range in 
each plot. For the r elevant radial ranges, we have also included the 
14 XMM groups of'Finog uenov et al T( l2006ll2007' ) that are not part 
of our sample and for which results have been extracted in a similar 
way. There is convincing evidence of an increasing Fe abundance 
with system temperature in the cores of these groups, but the trend 
progressively disappears as one moves to larger radii, and is sig- 
nificant at less than la beyond O.Srsoo- This conclusion is only 
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strengthened by the inclusion of the Finoguenov et al. data points 
in the correlation analyses, with which rises to +4.1 within 
O.lrsoo and drops to +2.1 within 0.1 r ^ O.Srsoo. The overall 
picture for Si is very similar, with a clear positive correlation within 
O.lrsoo (o"if = +2.5) but none further out. We emphasize that the 
absence of a clear trend at large radii supports our use of a sin- 
gle parametrization to describe the abundance profiles beyond the 
radius of X-ray detection where n ecessary. Outs i de the core, sim- 
ilar results have been reported by iTamura et al.l ( |2004|) for T and 
Z derived within the range 50 -200 kpc in a combi ned sample 
of groups and clusters, and bv lFinoguenov et aT lilOOl") for Z-pa of 
T ^ 1 keV systems derived at 0.2ri8o. However, at small radii the 
clear trend seen here differs from results for clusters; for example 
lOe Grandi et alj ( 12004) find no clear correlation between Z and T 
measured inside r25oo for their T > 3 keV systems. 

The observed trend in the group cores could potentially by 
induced by systematics; as suggested by the work of Buote ( 2000a), 
a Z-T correlation could arise due to the Fe and Si biases, which are 
more pronounced at low T. However, for most groups, the results 
in the core are based on mass-weighting spectral results extracted 
within narrow radial ranges in which T does not vary much. In the 
few cases where it does, two-temperature models were employed 
in spectral fitting, as described in Paper I, suggesting that we have 
generally resolved any temperature structure that could be affecting 
the abundance determinations. Further, the same overall trends are 
seen for the emission-weighted abundances, using Z directly from 
spectra extracted across the relevant radial range (which in some 
cases require two-temperature models for a satisfactory fit). This 
argues against a systematic effect related to the modelling of the 
density profiles in the group cores. We therefore conclude that the 
inferred trends should be reasonably robust against such biases. 



3.3 Supernova feedback 

In the context of elucidating the nature and history of galactic feed- 
back in the groups, the SN contribution to ICM energetics is an 
important quantity. Based on the adopted SN yields, we computed 
the energy _Esn released by the two SN types, assuming that energy 
has been released to the ICM in proportion to the metals at a level 
of 10^^ erg per supernova. Under this simplifying assumption, the 
specific SN energy - i.e. the total energy normalised by the num- 
ber of ICM particles - released within rsoo is shown as a function 
of grou p temperature in F ig. [TTI In a study of cool clusters out to 
0.4risn. rFinoguenov et al.l ( l200ll ) found a clear increase in specific 
SsN with T across the range T = 1-3 keV, with the trend flat- 
tening above T > 3 keV at 1.2 keV particle^ ^ (shown as a dotted 
line in Fig.ll U. Our results sample the low-T end better and reveal 
no significant trend within rsoo = +0.1), although a positive 
trend of weak significance {gk ~ +1) is found within O.Grsoo ~ 
0.4ri8o. Mainly as a consequence of the factor of two lower abun- 
dances in outskirts, the sample mean o f -Esn is correspondingly 
lower for our groups than found for the iFinoguenov et alj j200lh 
clusters. Note that Fig.[TT]also suggests that iJsN becomes compa- 
rable to the thermal energy iSth = (3/2) fcT of the ICM particles 
below r ~ 0.5 keV. While this is in excellent agreement with re- 
sults of cosmological simulations incorpo rating galactic outflows 
jPave. Op penheimer & Sivanandam' 2008), our results are incon- 
sistent with the specific Bsn declining with (T) as predicted by 
those simulations. 

The total SN energy released in the groups spans the range 4 x 
10)^^-2 X 10®^ erg. Although SN la contribute about ~ 20 per cent 
of the iJsN within the group core, the vast majority (~ 95 per cent) 
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Figure 11. Released SN energy per ICM particle inside rsoo- Dashed line 
marks our mean value of i? fa 0.6 keV p article"^ , while dotted li ne shows 
the mean found for T > 3 keV clusters jFinoguenov et al.ll200lh . Shaded 
area shows the region where Esti exceeds the thermal energy Bth of the 
ICM. 



of the total energy released within 7-500 has been provided by SN II. 
We emphasize however, that these numbers, as well as those shown 
in Fig. [TT] simply represent the SN explosion energy associated 
with the metals in the ICM and not necessarily the energy deliv- 
ered by SN to the hot gas. The two quantities would differ if, for 
example, the fraction of synthesized metals reaching the ICM is dif- 
ferent from the corresponding value for the SN explosion energy. 
We will return to these issues in Section l42l noting for now that 
uncertainties on the inferred i?sN are not necessarily dominated by 
those associated with our gas mass estimates or the choice of SN 
model yields. 



3.4 The origin of the central iron excess 

The presence of a central Fe excess in these groups parallels results 
for clusters w hich show this feature t o be common in cool -core 
systems (e.g., IFinoguenov et alj|200d IPe Grandi et al.ll2004l) . As 
for clusters, the excess can extend well beyond the optical extent of 
the central galaxy and can show a dip at the very centre, inside r < 
0.5-D25, where D25 is the ellipse outlining a B-band isophote of 25 
mag arcsec"^. As shown in Fig. [6] the enrichment in group cores 
has not been exclusively provided by SN la, but the Fe excess itself, 
relative to the Fe level at large radii, has largely been provided by 
SN la. This is broadly consistent with the idea that the central early- 
type galaxy is responsible for the Fe excess. Wh ile evidence for 
this conclusion is also obse rved in clusters (e.g. iBohringer et al.l 
|2004 iDe Grandi et aklliooj) . this scenario has not been tested in 
any detail in groups. One could a priori expect a similar situation 
to hold for lower-mass systems, since the fractional contribution of 
the brightest galaxy to th e total light is generally larger in these than 
in clusters ILin & Mohrll2004h . 

In order to investigate this, we computed the mass of Fe pro- 
duced by SN la within O.lrsoo. On average, this corresponds to 
a radius of r ~ 50 kpc, roughly matches the radial extent of the 
cool core where present (cf. Paper I), and corresponds to twice the 
maximal D2S extent of the central galaxy and to the radius where 
the radially weakly declining SN II contribution to Fe production 
becomes consistent with the level in the outskirts (Fig.[6]l. Results 
were normalised by the i^-band luminosity of the central galaxy 
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Figure 12. (a) Mass of Fe produced by SN la within O.lrsoo, normalised 
to tlie X -band luminosity of the central group galaxy. Shaded regions out- 
line the errors on the typical time-scales required for SN la to produce the 
observed Myc/Lk ratios, assuming a constant SN la rate with redshift. 
(b) The corresponding enrichment times, including also the effect of stellar 
mass loss in the BGG, and that of assuming an SN la rate evolving with red- 
shift according to equation @. Lines connect the time-scales for each group 
as computed under these different assumptions. A typical en'or bar (dom- 
inated by uncertainties in the local SN la rate) is shown at T = 0.3 keV. 
Dashed line mai'ks an enrichment time-scale of 14 Gyr. 



from Table [2] Assuming our adopted SN la yields along with the 
observed present-day SN la rate in nearby e arly-type (E/SO) galax - 
ies of i?o = 3.51J ? x 10"^" yr"^ L^^„ jMannucci et"aLll2005h . 



to = J-oii.i X iu - yr i^^ g 
the resulting time-scales required for SN la in the central galaxy to 
generate the observed Fe mass are illustrated in Fig. I12t t. Results 
from Si broadly agree with those from Fe. As can be seen, there 
is a clear dependence of the Fe M / L ratio, and hence required en- 
richment time-scale, on (T) within O.lrsoo, qualitatively mirroring 
the global behaviour inside rsoo (cf. Fig. [Sj. Both the mean and 
median enrichment time-scales for the sample are ~ 5 Gyr, which, 
given the increasing time-scale with T, is entirely consistent with 
the corresponding results of Bohringer et al. (2004), who find that 
times of > 5 Gyr are required in clusters. 

The immediate implication is that in many of these groups, 
SN la in the central group galaxy could well be responsible for 
the Fe excess seen within the cool core. Some caveats to this in- 
terpretation, however, include the facts that our method does not 
include the Fe that has become locked in stars over the relevant 
time-scales (although this might be a small fraction if star forma- 
tion largely ceased at high redshift in these central galaxies), nor do 
we account for any Fe expelled beyond O.lrsoo by galactic winds. 



or for the growth of stellar mass (~ Lk) of the BGG over the con- 
siderable time-scales involved. The inferred time-scales may thus 
be viewed as lower limits under the adopted assumptions. On the 
other hand, there could be a contribution to central enrichment due 
to stellar mass loss in the BGG and due to gas lost to group cores 
from satellites via galactic winds or environmental interactions. A 
more sophisticated approach would thus include the contribution 
from stellar winds along with a model of the stellar mass evolu- 
tion of the BGG. An argument for not considering the latter effect 
here, however, is the observation that the bulk of the stellar mass in 
many central clu ster ellipticals was in place already at z = 2 — 3 
(e.g., IStottet Zi r2008). so this may not play an important role. A 
case specific to our sample is NGC 1407, in which the BGG ap- 

S iears to have forme d over half its stellar mass already by z « 5 
Spolaor et al.ll2008l) . 

Mass loss due to stellar winds can be trivially incorporated in 
our calculations, and it is instructive to consider the impact of this. 
For this purpose, we consider all the Fe present within O.lrsoo (i.e. 
not just that produced by SN la), but subtract a uniform base level 
of 0.1 Z0 as seen in the group outskirts to obtai n an estimate of 
the excess iron mass within this radius. Following iBohringer et al.l 
([2004), the stellar mass loss rate M« from the central early-type is 
then estimated from 



M, 



1.5 X 10 



-11 Li 



(3) 



M0 yr ^ L© \ 15 Gyr 

where t, is the 'age' of the galaxy i iCiottietalJ[l99ll) . The rate 
of Fe loss is obtained by multiplying this estimate by the mean 
Fe mass fraction in stellar winds, which in turn depends on stellar 
metallicity. Assuming mildly supersolar metallicities for the BGG 
stellar population (Z = 1.5 Zp, corre sponding to [Z/H] ~ +0.2; 
see e.g. lSanchez-Bla zguez et al.ll2007 ) along with our adopted So- 
lar abundance table ( Greve sse & Sauval 19981) a nd a typical galaxy 
age of 10 Gyr ( Sanchez- Blazquez et alj 200^, the prediction of 
equation ([3} is added to the SN la iron production rate and the as- 
sociated enrichment times plotted in Fig.ll2b. 

On average, the inferred stellar mass loss can account for 
~ 25 per cent of the Fe within O.lrsoo, reducing the required en- 
richment times correspondingly. Nominally, this brings the esti- 
mated time-scales below 14 Gyr for all groups, alleviating - but not 
obviating - the need for additional sources of metal enrichment in 
the cores of high-T groups. Hence, subject to the adopted assump- 
tions, the time-scales may still be prohibitively long for systems 
with high M-pc/ Lk ratios in the core. Tentative evidence that satel- 
lite galaxies may have contributed to building up the central metal 
excess es in at least one of our groups is provided bv lMendel et al.l 
( l2009i) . who demonstrate a good correspondence between the ra- 
dial X-ray Fe distribution and that of stellar mass density in the 
NGC 5044 group out to r ~ 100 kpc, corresponding to five times 
the D25 extent of the BGG. Yet, barring a significant additional 
contribution to central Fe enrichment from other group galaxies, 
we are compelled to conclude that our approach generally under- 
estimates the metal production of the BGG over cosmic time. Pos- 
sible loopholes include invoking an SN la rate per unit Lk that 
evolves with redshift, similar to wh at has be en propos ed for early- 
types in more massive clusters (e.g.. lRenzin i 1993; Bo hringer et al.l 
l2004h . or a significant contribution to central ICM enrichment ei- 
ther from intracluster stars (see Sec. 14.21 ) or from th e sinking o f 
highly enriched low-entropy gas towards group cores jCorall2006h . 

We can briefly explore the former possibility, the effect of a 
redshift-dependent SN la rate. To this end, we assume an evolving 
SN la rate ria(f) per rest-frame time t and co-moving volume as 
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described by the star formation rate SFR(t) at a given time con- 
volved with the SN la delay-time distribution 4>(t), 

riaW (X I S¥K(t')^{t-t')At', (4) 
Jta 

where to is the age of the Universe at the initial redshift of star 
forma tion, here assumed to be at 2/ — 10. Following lEttoril 
(l2005h. SFRft ) is ad opted from the extinction-corrected model of 
IStrolger et ail j2004h , and <& is assumed to be a Gaussian of char- 
acteristic delay time r = 4 Gyr and g = 0.2 r (the result s are not 
highly sensitive to the exact choice of $; see lEttorill2005h . In this 
model, the SN la rate peaks at z « 0.8 at « 4.8 times the present- 
day value. Excluding here the effects of stellar winds (for straight- 
forward comparison to the results assuming a constant SN la rate 
with no stellar wind contribution), the resulting enrichment time- 
scales are plotted as empty squares in Fig. [T2b . This modification 
is seen to reduce the enrichment time-scales to ~ 5-7 Gyr even 
for the hot test groups, in agreem ent with similar estimates for clus- 
ter BCGs teohringer et al.ll2004h . Thus, in this model, the central 
Fe excess can easily have been generated by the BGG itself within 
half a Hubble time, without the need for additional sources of cen- 
tral enrichment. 

A final point worth mentioning here, related to mass loss 
through stellar winds, involves the associated abundance ratios. 
Since neither Si nor Fe is processed in the intermediate-mass stars 
primarily responsible for such winds, the a/Fe ratio of these winds 
(and of their parent stars) will be identical to that of the gas from 
which the stars formed. This ratio can generally be expected to be 
high in early-type galaxies, where star formation was largely ter- 
minated at high redshift as indicated by the supersolar a/Fe ratios 
typically found for the stellar population in these galaxies (e.g., 
ISanchez-Blazguez et alj|2007l) . Consequently, stellar winds from 
our BGGs should have a/Fe ratios approaching those of SN II 
ejecta. Interestingly, Fig. |6] indicates a mild central enhancement 
in SN II abundances which is at least consistent with a signifi- 
cant contribution from stellar mass loss associated with the BGG. 
Whether this feature can be entirely explained by stellar mass loss 
is a complicated question, however. A contribution from SN II to 
central enrichment is certainly still required, as the predicted sam- 
ple mean of Zsi/^Fc ~ 0.95 from SN la and stellar winds alone 
falls short of the observed mass-weighted mean of Zsi/Z^c ~ 1.10 
inside O.lrsoo. The radial extent of the SN II enhancement sug- 
gests either a contribution from stars not currently in the BGG, i.e. 
other group members or intracluster stars, or the action of some 
gas and metal mixing within ~ O.Srsoo. Moreover, the spatial 
distribution of stellar mass loss products is not straightforward to 
predict, since both [Z/H] and [a/Fe] in early-type stellar popula- 
tions tend to show radial gradients; while [Z/H] generally peaks to- 
wards the galactic centre, a/Fe gradients may be positive, negative, 
or fl at, and the mean absolute value may depend on galaxy mass 
jSan chez-Blazquez et al. 2007; Spolaor et al. 2009), implying that 
the Fe fraction and abundance ratio in stellar winds may be a highly 
non-trivial function of radius and of the environment of the BGG. 
A detailed exploration of this would likely require X-ray data of a 
superior quality and is beyond the scope of this study. 



4 ORIGIN AND FATE OF ICM METALS IN GROUPS 

We can briefly summarize the results of the previous Section as 
follows. Cool-core groups tend to show similar levels of central 
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Figure 13. Stellar-to-baryon mass ratio inside rsoo ■ The shaded region out- 
lines the 1-(T range of iW*/(M* + Micm) = O.ieto 03 covered by the 
clusters of lEttoriH2003h . Inset shows the corresponding result (same units) 
with the contribution from the BGG excluded. 

enrichment to those of more massive systems, but abundance gra- 
dients outside the core are stronger, leading to lower abundances 
at large radii. Group outskirts are heavily dominated by SN II en- 
richment, with the total energy per ICM baryon injected by SN 
somewhat lower than in more massive systems. At the same time, 
the global gas mass for a given amount of stellar light is lower in 
cooler systems within rsoo, leading to lower global metal mass-to- 
light ratios, both for SN la and SN II products. The latter result in 
particular may hold important information on the enrichment his- 
tory of these systems. In this Section we investigate a range of pos- 
sible explanations for this while acknowledging that this is not nec- 
essarily an exhaustive list. Fundamentally, these explanations can 
be divided into mechanisms causing (i) less efficient metal genera- 
tion or metal transfer to the ICM in low-mass systems, (ii) removal 
of metals from groups over cosmic time, or (iii) cosmological ac- 
cretion of less enriched gas in small systems. 

4.1 Galaxy formation efficiency and star formation history 

As an example of a mechanism belonging to the first category men- 
tioned above, we first consider the possibility of a reduced galaxy 
formation efficiency at low (T). Everything else being equal, this 
would lower the amount of Fe released per unit ICM mass, and al- 
though this caimot in itself explain the observed IMLR variations 
(unless the stellar initial mass function also varies with environ- 
ment), it could help explain the globally lower metal abundances 
seen in cooler systems. 

We can test this scenario on the basis of the gas mass-to- 
light ratio in Fig. |7] or equivalently, the ratio of stellar mass M» 
to that of baryons, i\fbar ~ -\- Micm (sometimes denoted the 
cold fraction /c, or the star formation efficiency) shown in Fig. 1131 
This plot indicates a general decrease in cold fraction with (T) 
{ffK = —1.5, strengthening to gk = —2.2 if excluding the some- 
what peculiar NGC2300; see Paper I), but with a tendency for 
this quantity to be higher than corresponding results for clusters 
(Ettori 2003). Hence, the galaxy formation efficiency actually ap- 
pears to rise in lower-mass systems within our sample, i n concor- 
dance with other observat i onal and theoretical finding s ("Lin et al] 
200i iRoineo et alj l2006l : iGonzalez. Zaritskv & Zabl udoff 200^; 
iDave et alj|2008l) . Hence, there is no relative dearth of stellar mass 
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within our groups, and we can readily exclude more exotic scenar- 
ios to account for the lower IMLR such as the possibility that the 
stellar population has simply produced systematically fewer metals 
in cool groups; the metallicity of the stellar population is typically 
at least ~ solar in many of the bright early-types in our groups, 
with no obviou s differenc es from central clust er galaxies (see, e.g., 
iHumphrev .^uote 2006; lBrough et al.ll2007b . However, it is also 
clear that the moderate trend of decreasing cold fraction with T is 
largely driven by the BGG (see inset in Fig.|13t. If considering the 
contribution to M,/Mbar from satellite galaxies alone, there is no 
systematic variation with group temperature (|o-_fs-j = 0.3, regard- 
less of whether NGC 2300 is included), and most groups are then 
consistent with the corresponding cluster results. 

For a fixed baryon fraction, a higher cold fraction in low-mass 
systems is perhaps no surprise. Cooler systems have lower char- 
acteristic entropies ('Ponman, Sanderson & Finoguenov' 200?) and 
hence cooling times (Sanderson, Ponman & O'SuUivan 2006), and 
so a larger fraction of t heir hot gas can be expected to cool out. 
As discussed by iGonzalez et alj j2007i) . a higher galaxy/star for- 
mation efficiency in low-mass systems could also be promoted by 
an increased efficiency of tidal interactions (boosting star formation 
rates) and a decreased efficiency of strangulation, i.e the stripping 
of any gas halo of late- type galaxies. 

It is nevertheless interesting to note that, if accounting also for 
the iron locked in stars, a higher fc in cool systems could poten- 
tially maintain a total IMLR comparable to that seen in clusters. 
However, it seems that only the stellar population of the BGG it- 
self could likely provide any significant modifications to a trend 
in this total IMLR with T, because Fig. [13] shows that the stellar 
metal contribution from satellite galaxies would not vary system- 
atically with T unless stellar metallicity does. For the sake of ar- 
gument, we note that any correlation of total IMLR with (T) can 
be made to vanish for our sample provided that the stellar pop- 
ulation in the BGG has mildly supersolar metallicity, [Fe/H] ~ 
+0. 2 in all g roups. Whil e not an unreasonable assumption (see, 
e.g., iHumphrev & Buot j l2006h . the resulting mean total IMLR 
of ~ 7 X 10"'^ would still be a factor of ~ 3 below the cor- 
responding value for clusters (w hich is ~ 0.02 for our adopted 
abundance table; iRenzinil |2004) . The contribution from satellite 
galaxies can raise this mean, but restores the temperature depen- 
dence, unless the (mass-weighted) mean stellar metallicity declines 
from Zfs ~ 2 Zq to ~ solar across our temperature range. In 
other words, a higher fc in cooler groups can only keep the to- 
tal IMLR independent of T and in accordance with cluster val- 
ues while maintaining the trends in Fig.[8]if a much larger fraction 
of synthesized metals remains locked in stars in cool systems. Al- 
though we lack the data to exclude this with confidence, it seems 
unlikely for several reasons, in part because it requires consider- 
able systematic variation in star formation histories and/or IMF 
and stellar binary fraction among the bright early-types in cool 
and hot groups. Secondly, a higher average stellar metallicity in 
cooler groups ~ in which the average galaxy stellar mass is lower, 
as discussed below - would conflict with the well-established local 
mass-metallic ity relation, which ap pears fairly insensitive to envi- 
ronment (e.g.. lBernardi et alj2003h . 

Including the contribution of any intracluster stars in our 
groups could help explain the trends in Fig. [T3] if such a com- 
ponent is more prom inent in poorer systems as suggested by 
[Gonzalez et alj ( |2007|) . However, given that its estimated contribu- 
tion to the stellar mass in clusters is still of order ~ 30 per cent 
(see Section [42), it seems improbable that it can account for the 
order-of-magnitude variation in i\fgas/L_g among our groups, nor 
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would it explain the seriously reduced IMLRs in cool systems. In 
summary, it seems clear that a systematically varying galaxy for- 
mation efficiency, or even a contribution from stars not included 
in our estimates of L_g or M,, cannot account for the results in 
Fig-ID Given the higher cold fraction in cooler systems, one would 
naively expect higher ICM abundances in these, at variance with 
our results. 

Although a varying star formation efficiency cannot explain 
the trends in Fig. [8] it is still instructive to also explore the po- 
tential role of variations in average star formation history among 
the groups. For example, is it conceivable that ICM enrichment 
has somehow been sufficiently 'delayed' in lower-mass systems 
to help explain their lower IMLR? To investigate this, we show 
in Fig. [14] the mass-weighted Si/Fe ratio of the groups as a func- 
tion of the corresponding Fe abundance, both restricted to radii 
r ^ O.Srsoo inside which we have reliable measurements for all 15 
groups. Systems have been grouped according to their mean tem- 
perature. Nominally, the data show a negative correlation with Zpc 
at the 2.0-a level, suggesting a higher Si/Fe and lower overall Fe 
abundance in cooler groups within these radii. The statistical errors 
are large, however, clearly exceeding the systematic variations, so 
the results have also been accumulated in three temperature bins 
for clarity, as illustrated in the figure inset. This confirms the initial 
impression, demonstrating that the ICM in cooler and less enriched 
groups tends to contain a relative preponderance of SN II ejecta. 

In analogy with stellar population studies, the tendency for the 
ICM to be less Fe-enriched and display higher a/Fe ratios in cooler 
groups suggests a chemically 'younger' ICM in these groups. On 
average, galaxies also tend to be slightly bluer in the cooler groups, 
as indicated in the figure inset, and furthermore are less massive, 
with the average stellar mass from equation l[2j showing a positive 
correlation with (T) at 3.0 and 3.4-a significance for BGGs and 
satellites, respectively. While several explanations may be possible 
for these results, the data are at least consistent with a scenario in 
which the (generally bluer and less massive) galaxies in our cooler 
groups have experienced slightly more extended star formation his- 
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tones, whereas the conversion of the available gas into stars was 
completed earlier in the more massive galaxies in hotter groups. As 
a consequence, star formation and metal production (from SN la 
in particular) might be progressively 'delayed' in lower-mass sys- 
tems. 

Although these results should be regarded as tentative {ok ~ 
— 1.3 for a correlation between Zsi/Zpa and total B-K colour), 
they could be a manifestation of the observation that more mas- 
sive early-type galaxies tend to be older, i.e. have completed the 
bulk of their star formation (and the associated ICM enrichment 
by SN II a nd stellar winds) at higher redshift ('downsizing', e.g., 
ICattaneo"e t al. 2008), as well as an effect by which such galaxies at 
fixe d mass tend to be slightly older in more massive environments 
(e.g. lSanchez-Blazguez et aTll2006i) . perhaps indicative of an envi- 
ronmental dependence of d ownsizing. This inte rpretation is also in 
accord with the inference of lLoewensteinI ( l2006h that the fraction of 
baryons processed through stars is lower in lower-density environ- 
ments compared to in rich clusters. Nevertheless, the inferred aver- 
age colour variations between galaxies in the cool and hot groups 
are small, A{B — K) < 0.3, so it remains highly questionable 
whether average star formation histories have varied sufficiently 
with group mass to explain the observed order-of-magnitude varia- 
tion in IMLR from both SN la and SN II. 

In conclusion, it seems unlikely that arguments invoking 
changing star formation efficiencies or histories with group mass 
can themselves explain the results in Fig. [8] The above results are, 
on the other hand, consistent with, for example, substantial gas (and 
hence metal) loss from within rspp in lower-mass systems. In fact, 
the simulations of lOave et alj ( |2008|) indicate that the higher cold 
fraction in cool systems (which, naively, could be expected to result 
in higher ICM metal abundances), is compensated for by galactic 
outflows that drive more metals out of the potentials of cooler sys- 
tems. We consider this possibility in Section l4!4l 

4.2 Metal release efficiency 

Metals can be released into the ICM through several mecha- 
nisms, including galactic outflows, stripping of enriched material 
by galaxy-galaxy and galaxy-ICM interactions, and direct in situ 
enrichment by intracluster stars. Here we consider these mecha- 
nisms in turn, posing the question of whether metals have been re- 
leased less efficiently into the ICM in cooler groups. One immedi- 
ately attractive feature of this possibility is its potential to account 
for the lower metal mass-to-light ratios of both SN la and SN II 
products. For example, gaseous stripping would likely affect both 
SN la and SN II products, and tidally stripped intracluster stars 
should release both SN la and SN II metals with equal efficiency. 
In addition, any metal release mechanism whose efficiency scales 
with system temperature might alleviate the need for a redshift- 
dependent SN la rate in the BOG (Section [3.4t . and help explain 
the observed Z-T trend in group cores. 



4.2.1 Galactic outflows 

Some fraction of all metals released by SN in the group members 
will never mix with the ICM but will remain in the parent galaxy 
and potentially participate in subsequent star formation. Fig.[8]then 
raises the possibility that SN ejecta were released less efficiently 
into the ICM in cooler systems. To evaluate the feasibility of this 
scenario, we can compare the total SN energy inferred from the 
ICM abundance pattern £sn,icm (Fig. lilt to that expected from 
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Figure 15. Ratio of SN energy released in the groups as inferred from 
ICM abundances to that expected from the stellar mass of the group mem- 
bers. Sh aded regio n mar ks the estimate for the T Ri 2-6 keV clusters of 
ISivanandam et alj feoogi) . 

the stellar population within a group, iJsN,*- Naively, one would 
expect rj — Ssn.icm/Ssn,* ^ 1, the exact value depending on 
the mass fraction of SN ejecta retained by the galaxies. 

For each group, Esn,* can be estimated from the stellar mass 
of the group members, noting that a Salpeter IMF with an as- 
sumed minimum mass for core-collapse progenitors of 8 Mq im- 
plies one core-collapse SN per ~ 150 M© of stars formed. In ad- 
dition, a typical stellar afPe ratio of ~ solar in the group mem- 
bers (Humphrev & Buote 2006) implies that at least ^ 0.3 SN la 
must have exploded for each SN II. Adding these numbers yields 
£sN,» ~ 9 X 1O''^(M,/M0) erg , and the resulting comparison to 
i5sN,iCM within r^oo is plotted in Fig.[T5] With a resulting sample 
mean and l-o error of 77 = 0.68 ± 0.10, our numbers are in broad 
agreement with existing theoretical and observational estimates of 
the mass l oss fraction of galaxies, w hich range from ry ~ 0.4 — 0.8 
(see, e.g.. lSivanandam et al.ll2009l and references therein). In par- 
ticular, as illustrated in Fig. [15] the latter authors find a value of 
7? = 0.84^^014 within r^oo in clusters using an approach based on 
the observed SN la rate in cluster early-types, while noting that the 
fraction could be as low as 35 per cent if the assumed SN rates were 
increased to within their 2-o upper limit. 

Interestingly, while this suggests that ICM enrichment has 
been fairly efficient for our sample as a whole, there is also a 
tendency (ok ~ +1.7, rising to +2.5 if excluding the outlying 
NGC 2300) towards higher rj in hotter groups. This result could 
be anticipated from Figs.[TT]and[T3] Since the estimated Ssn.icm 
per ICM particle does not vary systematically with T, the trend in 
Fig.[T5]is a consequence of the systematic decrease in cold fraction 
with T. The result is consistent with the notion that SN products 
have been released more efficiently into the ICM in hotter systems, 
although other explanations are possible. These include the possi- 
bility of the stellar IMF and/or mean metallicity varying system- 
atically with T, as this would lead to systematic variations in the 
predicted number of SN per stellar mass and in the estimated stel- 
lar mass itself. 

Nominal values of rj exceeding unity are seen in a few cases, 
but all results are consistent with 77 < 1 even without accounting 
for uncertainties in SN rates and yields. Other possible effects in- 
clude the presence of optically undetected intracluster stars, which 
would contribute to i?sN,iCM only, and more speculatively the pos- 
sibility of pre-enrichment, whereby some of the ICM material ac- 
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creted by the groups was enriched by galaxies not currently in the 
group. In addition, we cannot exclude that ICM masses may be 
slightly overestimated in some groups, nor that the SN la contribu- 
tion to -BsN,* may be underestimated (cf. Sivanandam et al. 2009). 
The fact that r] is not constant and gene rally different fro m unity 
is no surprise though. As pointed out by Dave et ^ j2008l) . X-ray 
emission-weighted a/Fe ratios, used here to infer -Bsn.icm, may 
not be a good indicator of the relative number of SN la and SN II 
that have exploded within the group, but rather reflects the metal 
distribution mechanisms at work. 

A possible interpretation of the systematic trend in Fig.lISlis 
that SN-driven outflows have been more efficient at releasing met- 
als into the ICM in hotter groups. Whether this might be expected 
on the basis of the higher average galaxy mass in such groups is not 
immediately obvious though. A deeper galactic gravitational po- 
tential might disfavour the development of a freely outflowing star- 
burst wind, but, conversely, star formation may generally have pro- 
ceeded more rapidly and vigorously in massive galaxies. An alter- 
native explanation invokes the possibility of AGN-driven outflows 
helping to release metals into the ICM. Observational evidence for 
such outflows to exhibit significantly supersolar abundances for a 
rang e of elements has been reported ( Ar av et a l. 2007; Wang et al. 
12007) . and numerical simulations jMolf et alJi2007l) indicate that 
AGN outflows could contribute significantly to cluster enrichment. 

Since most of our groups show some evidence of ongoing 
or recent AGN activity (Paper I), the relevant question is whether 
such activity, and its impact on the ICM, scales with (T) as re- 
quired by Fig. [15] More specifically, since AGN activity may ex- 
pect to scale with stellar bulge mass of the host galaxy, we can test 
whether the trend in Fig. [15] is explicable on the basis of higher 
BGG bulge masses in hotter groups. Adopting again the distinc- 
tion between cool and hot groups introduced with Fig. [9] we note 
that 1] is significantly higher in hotter (T > 1.1 keV) groups by 
58± 18 per cent. Using the BGG central stellar velocity dispersions 
(To,, listed in Table [2] we find some tendency for these to increase 
systematically with T at the same level of significance as that of 
rj vs. T (aK = +1-7), However the trend is considerably flatter 
than that of r?(r), with results showing a mean and 1-cr error of 
((TO,.) = 267 ± 15 and 278 ± 10 km s~^ for cool and hot groups, 
respectively, indistinguishable at the l-cr level. Qualitatively similar 
results are obtained for {(Tq which is expec ted to scale with the 
mass of the central superm assive black hole jperrarese & Merrit^ 
l200d iGebhardt et alj|200oh . In other words, there are no indica- 
tions that central black hole masses are generally higher in hotter 
groups, and this result still ignores the fact that relatively more gas 
needs to be enriched in more massive systems. Hence, unless more 
massive black holes are considerably more efficient /or their mass 
at expelling metals from their host bulges, these considerations sug- 
gest that the trend in Fig.[T5]is not driven by increased AGN activity 
in hotter groups, at least not if associated with the BGG only. 

4.2.2 Stripping and interactions 

Gas and metal loss from galaxies induced by interactions with 
the ICM could also help explain the trend in Fig. [15] An indica- 
tion that some (presumably enriched) gas has been removed from 
our group galaxies by other means than SN outflows is provided 
by the factor of 2-3 deficiency in Hi mass derived for the spi- 
rals in NGC 5044 and NGC 7619 c ompared to field spirals of sim- 
ilar o ptical morphology and size jSengupta & BalasubramanvamI 
l2006h . In addition, evidence for ongoing or past stripping as- 
sociated with galaxy-ICM interactions is present in a few of 



our groups ( iRasmussen. Ponman & Mulchaevll2006l : lMendel et al] 
20 0S)- General support for this interpretation comes from the work 
of iFinoguenov et all j2006l) who find a larger scatter in metallic- 
ity than in entropy in the group cores within their sample, which 
they interpret as evidence for diversity in the origin of gas in the 
cores. On the other hand, iRenzinil ( l2004h notes that if, e.g., ram 
pressure stripping were important for ICM enrichment, then hotter 
clusters with higher galaxy velocity dispersion should show rela- 
tively higher ICM abu ndances. While this is no t observed in his 
results, nor in those of iBaumgartner et al.l ( l2005h . such a trend is 
seen for our sample, and, interestingly, is confined to the regions of 
highest ICM density (Fig.[TOl(. 

This trend is certainly consistent with a density-dependent 
metal release mechanism in group cores, whose efficiency scales 
with galaxy velocity dispersion. If ram pressure stripping is the 
culprit, we might expect the ICM abundance in the central re- 
gion to scale more strongly with the characteristic ram pressure 
in this region than with group temperature. In order to test this pos- 
sibility, the velocity dispersion of each group (Table [TJ, was 
used to estimate the characteristic ram pressure within O.lrsoo- 
Assuming P^am = (Mgas/V)(T^, with V = (4/3)7^(0.1r■5oo)^ 
we find evidence for a moderately significant correlation between 
Pram and the mass-weighted Fe abundance within O.lrsoo ((tk — 
-1-2.1), though the correlation of Zp^ with (T) is actually stronger 
[aK = +2.6 for own grou ps, -1-3.0 when including those of 
IFinoguenov et alj l l2006l l2007^ as discussed in Section [T2l . In ad- 
dition, the corresponding correlation of Zsi with Pram is weak 
((TK — +1.1 compared to +2.5 for Zsi vs. (T)). While keeping 
in mind that the estimated may not always be representative of 
the typical 3-D galaxy velocity in each group if affected by velocity 
anisotropics or small-number statistics (although 13 out of our 15 
groups have more than 10 members), we thus find that results are 
consistent with ram pressure stripping contributing to enrichment in 
group core, but the evidence is tentative at best, and we cannot rule 
out other density-dependent mechanisms as more important. More- 
over, as this process should be inefficient at large radii (and where, 
indeed, no clear trend in Z with T exists), we conclude that the 
reduced efficiency of stripping through galaxy-ICM interactions in 
cooler systems is unlikely to account for their lower IMLRs or for 
the trend in Fig. 1151 

Note that if the metallicity-temperature trend in group cores 
does reflects an increased efficiency of environment-driven gas loss 
from galaxies in hotter group s, then the absence of a similar trend 
in systems with T > 3 keV jPe Grandi et al.ll2004l) may suggest 
that such processes 'saturate' at T ~ 2-3 keV, perhaps as a conse- 
quence of near-complete gas removal from core galaxies in clusters 
above these temperatures. It has yet to be established to what ex- 
tent this possibility is in accord with simulation results, which sug- 
gest that only a minor fraction (~ 10 per cent) of the global ICM 
metal budget even in massive clusters can be accounted for by ram- 
pressure strip ping, though with the fraction considerably higher in 
cluster cores jPomainko et al .120061 : iKapferer et al.ll2007h . 

Other interaction processes involve galaxy-galaxy interac- 
tions such as harassment, tidal interactions, and mergers. The latter 
process in particular may also release metals to the ICM, in part be- 
cause a significant fraction of a galaxy's ISM may shock-heat and 
become unbound during major mergers, and in part because the 
merger activit y may further tr i gger a starburst- or AGN-induced 
outflow (e.g.. Icox etai]|2006l : [Hopkins et al.ll2006h . Our data do 
not allow useful constraints on the contribution of this process to 
ICM enrichment, however, and cosmological simulations may still 
lack the dynamic range to robustly assess whether it has been sub- 
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stantially less efficient or influential in low-mass systems compared 
to clusters. Naively, one would expect tidal stripping and mergers 
to be more effective in poorer systems owing to the lower velocity 
dispersion, suggesting that these mechanisms cannot explain the 
trends in Figs.[8landll5l 

4.2.3 Enrichment by intracluster stars 

Intracluster (IC) stars represent a non-negligible component of 
the baryons in groups and clusters, and could pollute the ICM 
very efficiently since their ejecta are released directly into the 
ICM. Observations suggest that ~ 30 per cent of the stel- 
lar mass in clusters within rsoo may be contained in IC stars 
dOonzalez etai]|2007h . but it remains debatable whether this frac- 
tion de creases with total cluster mass as claimed bv lOonzalez et al.l 
( l2007h . or rises asymptotically with mass to a level of ~ 
50 per cent of M, as suggested b y other observational and the- 
oretical studies ( Lin & Mohr 2 0041: Conrov. Wechsler & Kravtsovl 
bOOllPurceilBuUock & Zentne]|2007h . 

Since this component is not accounted for in the cluster com- 
parison results in Figs. [7] and [8] including its contribution to La 
could aifect the apparent group-cluster discrepancies. For example, 
if IC light is less prominent in cooler systems as suggested by some 
of the above results, this could help explain the lower ICM abun- 
dances in groups. However, even a 50 per cent increase in IC stellar 
light in massive clusters is insufficient to bring the lower group 
IMLR within the factor of two of rich clusters that is explicable on 
the b asis of the lower large-ra dius abundances in groups. Further- 
more, [siv^Md^letal] j2009l) conclude that 20-25 per cent of the 
Fe in the ICM within rsoo can be accounted for by IC stars, im- 
plying that IC stars contribute in roughly equal proportions to the 
stellar mass and ICM Fe enrichment in clusters. This suggests that 
the resulting IMLR should not be significantly affected by the pres- 
ence of this component. Consequently, it seems we can exclude the 
possibility that intracluster light provides an important contribution 
to the trends in Fig. [8] This would be particularly true if this com- 
po nent is actually more important in cooler systems, as suggested 
bv lGonzalez et alj feOOvh . 

IC stars could still play an important role for ICM enrich- 
ment in group cores, however. A relevant question is whether the 
metal contribution from these can alleviate the need for a redshift- 
dependent SN la rate in the BGG to account for the Fe within 
the cool core (cf. the discussion of Fig. I12t . Where IC light in 
clusters can be separated from that of the central galaxy, the for- 
mer typically contributes ~ 80 per cent of their combined op- 
tical output within rsoo jSivanandam et alj |2009|) . This fraction 
will be lower inside O.lrsoo, because IC light is more widely dis- 
tributed than that of the BGG. For simplicity, if assuming that the 
IC light traces the group gravitational potential and that its K- 
ban d luminosity density is distribut ed according to an NFW pro- 
file n avarro, Frenk & Whitelll997h with a typical concentration 
c = 5-10, then 5-10 per cent of this light would fall inside O.lrsoo. 
The implication is that this component can potentially reduce the 
enrichment times in Fig. [12] by a factor ~ 1.4 (optimistically as- 
suming a total A'-band contribution of four times that of the BGG, 
of which 10 per cent falls within O.lrsoo). Although this factor as- 
sumes that none of this light is already accounted for in the adopted 
values of Lk.bgg, it may still be large enough to eliminate the 
need for an evolving BGG SN la rate, especially if the IC contribu- 
tion does increase with group mass. What IC stars probably cannot 
do, however, is account for the observed Z-T trend in group cores 
(Fig.llOt. Even if the fractional Fe contribution from IC stars within 



O.lrsoo increases by 50 per cent acros s our temperature range (at 
variance with the lGonzalez et al.1l2007l results), this would still fall 
an order of magnitude short of explaining the factor ~ 10 rise in 
Zpo within this radius. Arguments invoking a far superior enrich- 
ment efficiency from this component (compared to stars in galax- 
ies) would also fail, since Fig.[T5]suggests that galaxies are already 
fairly efficient in polluting the ICM. 



4.3 Cooling of enriched gas in group cores 

Turning now to the possibility that metals were removed from 
groups rather than not produced or released into the ICM, we first 
consider a scenario in which radiative cooling has removed en- 
riched gas from the X-ray phase. At a given fracti on of rsoo, cool- 

ing times tend to be lower in cooler systems (e.g.. lSanderson et al.l 
I — 

2006), so this explanation would be consistent with such systems 
exhibiting lower A/gas /is and missing relatively more metals, as 
suggested by Figs.|7]and[8] The process is nevertheless still limited 
by the observation of the central abundance excess, because cooling 
rates above a certain threshold would eventually destroy the central 
excess due to the associated inflow of less enriched material from 
larger radii. 

Radiative cooling would primarily deplete SN la products, due 
to their higher concentration at small radii where cooling is most 
rapid. Thus, the two questions we seek to address is whether the 
amount of Fe generated by SN la that is 'missing' in Fig.[8ji (or its 
A'-band analogue) can be accounted for by cooling, and whether 
that would ultimately destroy the central Fe excess in conflict with 
our observations. Since Figs. [8] and [T3] indicate that the properties 
of our most massive group, the T ~ 2.1 keV system NGC6338, 
are fairly representative of those of more massive clusters, we can 
use this group as a benchmark for characterizing deviations in cool 
systems from cluster results. On the basis of the A'-band IMLRia 
in NGC6338 deduced from Table[T] we evaluated the missing Fe 
mass in each group and hence the corresponding total missing gas 
mass for a given assumed Fe abundance of the cooled-out gas. Mak- 
ing the assumptions that the gas cooling out had Z ~ Z© as ob- 
served in the group cores, and, conservatively, that this gas has been 
cooling for 10 Gyr (corresponding to an initial redshift of z = 2 
for the adopted cosmology), resulting mean cooling rates could be 
computed. 

Fig. |16b shows the required cooling rates as a function of 
mean group temperature, demonstrating that modest average rates 
of a few tens of M© yr^^ would be necessary under the above 
assumptions. Note that NGC6338 itself has not been included in 
the plot, nor has the (T) ^ 1 keV system NGC4325, which ac- 
tually exhibits a slightly higher A-band IMLRia than NGC6338 
and so is not Fe-deficient according to our assumptions. Of course, 
the derived rates assume that gas in group cores was already en- 
riched to solar levels by z ~ 2 and has been subject to unimpeded 
cooling ever since. These estimates are thus likely representing a 
best-case scenario, and the total missing gas mass becomes corre- 
spondingly larger if assuming a lower Z^c for the removed gas. 
For example, adopting instead the current Mpe /A/gas ratio within 
rsoo in each group as representative of the assumed Fe content 
of the cooled-out gas, the required M rises to a sample mean of 
~ 130 M0 yr~^. Such values seem prohibitively large at present, 
but we cannot a priori exclude cooling rates of this magnitude in 
the past. Note also that cooling should be more efficient in lower-T 
systems, and although Fig.|16b indicates a weakly significant ten- 
dency {gk = +1.2) towards higher cooling rate requirements in 
hotter systems, there is a clear negative trend with T if normalising 
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Figure 16. The time-averaged gas cooling rate M required to account for 
the shortfall of IMLRia compared to NGC 6338, assuming cooling over a 
10 Gyr time-scale. Results are plotted as a function of (a) (T), (b) the B-K 
colour of the brightest group galaxy, (c) the star formation rate of the latter, 
where available (see text for details), and (d) the stellar mass of the latter. 
Dashed line in (c) represents equality, while lines in (d) mark out constant 
fractions of cooled-out gas mass relative to the stellar mass of the BOG, 
with fractions as labelled. Uncertainties on M reflect the relative en'ors on 
IMLRia. 



the required M by the total gas mass (cta- 



-2.6) or gravitating 



mass {uk = —2.1, assuming A/grav oc T^/^) of the group. 

Although the central temperature is not observed to drop be- 
low ~ 40 per cent of the peak value in any of the groups (Pa- 
per I), it is still relevant to ask whether there are any indirect signs 
of ongoing strong cooling where this is most needed. If central 
cooling is currently proceeding at the required rates, and some of 
the cooled-out gas participates in star formation, one might ex- 
pect systematically bluer BGGs in higher-A/ groups. We inves- 
tigate this in Fig. 1 16b . which reveals no significant dependence 
between M and BGG colour (aK ~ —0.7), not even if normal- 
ising required AI by Lk of the BGG (aK = —0.5). For the galax- 
ies with the relevant available data, we can be more specific and 
compare the required M to the current central star formation rate 
(SFR), noting that some low-level 'residual' star formation (SF) 
seems to be taking place in man y nearby early-type galaxies (e.g., 
ICombes. Young & BureajboOT ^. Hence, unless cooling and SF is 
highly episodic, one might again expect a correlation between M 
and current SFR. We have used available Ha and UV data from 
the literature to test this scenario (again excluding NGC 4325 and 
6338) . Ha-based SFRs from SDSS-DR4 data (Bri nchmann et al.l 
l2004h were adopted for NGC 4073 (the BGG in MKW4) and 
NGC 5846 , while Ha fluxes were ta ken from N ED fo r NGC 383 
and 741 jCapetti. Verdoes Kleiin & Chiaberg j |2005|), and UV 
fluxes from NED fo r NGC 1407 ( Rifatto. Longo & Capacciolil 
[l995) and NGC 4125 jDale et al. 2007). For th e latter foii r galax - 
ies, SFRs were estimated using the relations of iKermicutj jl998l) , 
which should be sufficiently accurate for our purposes. The UV 
luminosities have been corrected for Galactic extinction but not in- 
ternal dust extinction, which should be small for these early-types. 



The Hof fluxes typically derive from the central galaxy regions, so 
some AGN contribution is possible in these cases. For example, the 
Hq emission from NGC 4073 cannot be unambiguously classified 
as either due to star formation or an AGN, so we have treated this 
result as an upper limit. The results are plotted in Fig.ll6fc and show, 
if anything, a trend that runs opposite to expectations in our cool- 
ing scenario, and with SFRs generally several orders of magnitude 
below the required cooling rates. 

Thus, we find no convincing evidence that the required central 
cooling in the groups is associated with any recent SF in the BGG, 
with the necessary cooling rates generally higher than even current 
SFRs in central galaxies of 'cooling flow' clusters (t hese rarely ex- 
ceed 5-10 Mq yr~^ and are often much smaller; cf. lRaffertv et al.l 
l2006h . This is corroborated by the fact that there is no clear ten- 
dency for the BGG to be bluer in cooler systems (aK ~ —0.4) 
where cooling is expected to be more efficient. Further, large 
amounts of recent cooling-induced star formation seems incompat- 
ible with the inferred ages of the stella r population in many of th e 
BGGs, which are of order 10 Gyr (e.g.. lHumphrev & Buotdl2006h . 
This issue can potentially be circumvented in the - perhaps some- 
what unlikely - event that much of the cooled-out gas has been de- 
posited without participating in star formation. If so there should be 
substantial (several 10^*^ M©) reservoirs of cold gas residing in the 
group cores, but Hi measurements listed in NED, where available, 
only reveal upper limits, with NGC 5846 as the only exception. 

The above considerations argue against significant recent 
cooling, and reveal no indications that ongoing cooling is more im- 
portant where required. However, this does not preclude the possi- 
bility of stronger cooling activity in the past, where it could have 
contributed to building the current stellar mass of the BGG. If so, 
and if cooling is indeed responsible for removing metals, one might 
expect the stellar mass Af,,BGG of the central galaxy, and the ra- 
tio of this to total ICM mass within rsoo, to correlate positively 
with the required amount of cooled-out gas. This is not clearly the 
case though, with correlation significances of only ok = +1.3 
and +0.1, respectively. An example is demonstrated in Fig. |16ti, 
which compares AI to Af,,BGG. In a few cases, the required cooled 
gas mass nominally represents ~ 100 per cent of the current stel- 
lar mass of the BGG (implying that the vast majority of the metals 
proposed to have cooled out cannot have been supplied by the BGG 
itself), while on average it corresponds to about half this value. 

Although the latter case may not necessarily be completely 
at variance with the typical history of mass build-up and star for- 
mation in these galaxies (which may still have assembled a signif- 
icant fraction of their current stellar mass since z ~ 2), there is 
an additional and very powerful constraint on the allowing cool- 
ing rates arising from the presence of the central Fe excesses. As 
the typical estimated iron injection rate from the BGG is of or- 
der ^ 0.02 Mq yr~^ if including both SN la and stellar winds, 
the cooling rates implied by Fig. [16] would destroy the observed 
excesses on typical time-scales of only ~ 10^ yr through the as- 
sociated inflow of less enriched material from larger radii (while 
also flattening the central entropy distribution for the same reasons, 
although a discussion of this is beyond the scope of the present 
study). This severely compromises the viability of recent cooling 
as a mechanism to reduce the global X-ray iron mass-to-light ratios 
in cooler systems. 

Could central cooling and inflow of low-Z gas still play a role 
for establishing the central Z-T trend in Fig. HOP If repeating the 
cooling calculations underlying Fig. [T6] using the 'missing' metal 
mass inside O.lrsoo (again assuming that IMLRia in NGC 6338 is 
representative of the case of no missing metals), then the typical re- 
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quired cooling rates are a factor of ~ 3 smaller, with some tendency 
for the required M/Mgaa to decrease with T {<jk = —2.5) and for 
M/Lk.bgg (fTif = —1.7) to decrease with (_B-A')bgg- This sug- 
gests that central cooling could have played some role in removing 
metals from group cores and suppress the central Fe excess in the 
coolest systems, although, again, unimpeded cooling at the rates 
required to explain the Z-T trend within O.lrsoo in Fig.llOlwould 
generally have destroyed the observed Fe excesses. 

In summary, although cooling rates under optimistic assump- 
tions may be sufficiently low in some systems for this mechanism 
to remain a viable means of removing some metals from the X-ray 
phase, the rates must generally be far below those required to ac- 
count for the global shortfall of Fe within rsoo implied by Fig.[8l In 
addition, there are no clear indications that cooling is, or has been, 
more prominent in the systems with the largest deficit of SN la 
metals. As mentioned, cooling would also primarily affect SN la 
products, thus likely still requiring a separate mechanism to also 
account for the lower IMLRn in cooler systems. We thus conclude 
that central cooling cannot generally have been important in estab- 
lishing the low IMLR of groups. 

4.4 Recent gas and metal loss from groups 

Another way to induce a lower IMLR in groups is through gas or 
metal loss from the potential wells of low-mass systems. This could 
help account for the lower gas mass-to-light ratios (via bulk gas 
loss) and potentially also the lower abundances (via preferential 
loss of highly enriched gas). As a means of discriminating between 
bulk gas loss and selective loss of enriched gas, we note that only 
the latter process would generally reduce the inferred abundances. 
The efficiency of such a mechanism should depend on group tem- 
perature, and our sample is large enough to allow a crude test for 
such an effect as demonstrated in Figs.|9land|10| 

Beyond O.Srsoo, Fe has mainly been produced by SN II 
(Fig. |6l(, so Zfc should here provide a crude measure of the abil- 
ity of the systems to retain the products of starburst-driven galactic 
winds. There is little evidence from Fig.[TO]that this ability varies 
systematically across the temperature range probed here {gk < 1), 
although the correlation significance rises to gk ~ +1.2 if exclud- 
ing the fossil group NGC741 which is extremely iron-deficient in 
the outskirts. Hence, although statistics at large radii are limited 
and dominated by the hotter groups in the sample, with only three 
systems from the 'cool' (T ^ 1.1 keV) subsample included in the 
last panel of Fig. (TO] it seems clear that selective removal of en- 
riched gas from cooler systems does not contribute significantly to 
the trends in Fig. [8] Moreover, given the strong radial decline of 
the mean Fe profile for the sample, a possible alternative - bulk gas 
loss from outskirts without preferential loss of metals - would pre- 
dominantly have removed low-Z gas, thus leaving material with a 
higher abundance behind. Hence, while such a mechanism could 
have lowered Mgaa/is in cool systems, it cannot account for the 
seriously reduced IMLRs in these environments. 

Gas and metals can certainly still have been lost from groups. 
The absence of a clear temperature-abundance trend within 0.5- 
1 -0 ''500 could potentially be an artefact of small-number statistics 
and does not preclude the possibility that enriched gas has been 
expelled well beyond 7-500 in all systems. Any such activity might 
be ascribed to galactic outflows driven by powerful SN or AGN 
activity, displacing ICM metals to radii where they would gener- 
ally escape detection in these low-T systems. In the absence of 
significant mixing, highly enriched high-entropy material released 
by such outflows would rise buoyantly in the ICM until in entropy 



equilibrium with the surroundings. We here consider the feasibility 
of such a scenario for both SN- and AGN-driven outflows. 

4.4.1 Supernova-driven winds 

Supernova-driven galactic winds should generally be dominated by 
SN II products. Indeed, observations of extraplanar X-ray emis- 
sion associated with hot gas around star-forming galaxies indicate 
significantly enhanced a/Fe ratios of 2-3 times the Solar value 
jStrickland et alj|2004h . The radial increase in the Si/Fe ratio seen 
for our groups outside the central galaxy is, in fact, only repro- 
duced in simulations where the feedback scheme favours the ejec- 
tion of SN II products from galaxies o ver those of SN la (see 
iTomatore et al.ll2004 iRomeo et alj|2006l) . This raises the interest- 
ing possibility that the radial increase in the Si/Fe ratio inferred for 
many of our groups may be explained by SN Il-dominated mate- 
rial that has risen in the ICM to mix with it at radii corresponding 
to a significant fraction of rsoo, and occasionally beyond. This rise 
could have occurred by means of buoyancy forces and/or bulk ki- 
netic energy impart ed to the gas. Since s tar-forming spirals tend to 
avoid group cores ( IHelsdon & Ponma n 2003), these galaxies will 
typically already be surrounded by high-entropy gas, so convec- 
tive displacement of the enriched ejecta to large radii may be in- 
efficient and mechanical work energetically favoured for removing 
large amounts of gas from within rsoo • 

Stringent tests of this possibility are impeded by the fact that 
SN wind properties remain poorly established. In particular, de- 
tailed hydrodynamical simulations indicate that much of the wind 
fluid in an SN-driven outflow may escape the galaxy in a high- 
entropy phase that would be difficult to detect with current X-ray 
instrumentation tStrickland & Stevens 2000). To obtain a crude es- 
timate of wind properties, we may consider the tail gas being lost 
from the NGC 2276 starburst galaxy, a me mber of the NGC 2300 
group in our sample (Rasm ussen et alj2006l) . This material is likely 
to contain a significant contribution from starburst outflows, though 
possibly with some hot ISM and surrounding ICM mixed in. The 
entropy of this gas is ~ 300 keV cm^, although this figure is proba- 
bly uncertain by a factor of ~ 2 due to the unknown geometry of the 
tail along the line of sight. Taken at face value, this gas would rise 
due to convection to at least r ~ 0.6rsoo if assuming the adopted 
ICM density profile for the group. This consideration suggests that 
starburst wind material can rise to significant radii in the ICM by 
means of buoyancy, provided that mixing with the surrounding ma- 
terial is suppressed at least initially, e.g., by magnetic fields. 

If SN Il-enriched material from starburst winds has been pref- 
erentially lost from lower-mass systems to generate some of the 
trends in Fig. [8] one would naively expect a positive correlation 
between Msi/M-pc and (T). While Fig. [14] already suggests oth- 
erwise, we examine this possibility more directly in Fig. [TT] again 
showing this ratio within O.Srsoo only, to avoid reliance on our 
temperature-independent abundance parametriz ations. The results 
are compared to the analogous values derived bv lSatoetalJj2007h 
for the T ^ 3-3.5 keV clusters AWM 7 and A1060 inside 0.35 
and 0.25rigo (« 0.6 and 0.4rsoo), respectively. Fig. 1 17| indicates 
that Si has not been preferentially lost from the cores of lower- 
mass systems for the T-range probed here. To the contrary, there 
is mild evidence for a negative correlation {gk — —1.8) although 
this is driven by the three coolest systems, without which any cor- 
relation vanishes (gk ~ —0.3). There is also no significant trend 
at larger radii (gk = +0.5), as evidenced by the inset showing 
the corresponding results within 0.5-0.7 rsoo for the nine groups 
with data extending across this range. Everything else being equal. 
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Figure 17. Ratio of silicon-to-iron mass inside O.SrsoQ. Also shown are 
the 90 per cent confidence intervals derived for the clusters AW M 7 (darker 
shade, dashed lines) and A 1060 (lighter shade, dotted lines) bv lSato et alj 
Inset shows the con'esponding values within the range 0.5-0.7 rsoo 
for the groups with data across this range. 



buoyancy-driven bubbles of SN II material should be able to rise to 
larger radii in cooler (lower-entropy) systems, but there is no evi- 
dence that the AIsi/M-pc ratio is higher at low (T) within this radial 
range. In other words, there is no clear evidence for a loss of SN II 
products from the central regions of cooler groups, nor for any such 
material to have been d isplaced to larger radii in those systems. 

Comparison to the lSato et al. 112003) results thus suggests that 
any selective removal of SN II products from the central regions of 
low-mass systems by starburst winds must have occurred with com- 
parable efficiency in systems with T < 3.5 keV, and so are unlikely 
to explain the trends seen for SN II products in Fig. [8] This conclu- 
sion is consisten t with the results s een fo r the combined group- 
cluster sample o f iFinoguenov et alj l l2000l) . which admittedly con- 
tains only four T < 2 keV systems. These authors find a roughly 
constant Zsi/Zpa ratio at r = 0.2ri8o for systems with a corre- 
sponding Si M/L ratio or gas mass fraction /gas below a certain 
threshold. Above this value, the Si/Fe ratio appears to rise dramati- 
cally, with the transition occurring around Si Ad/ L ~ 0.01 M^/Lp) 
or /g as ~ 0.1, corresponding to T ~3-^ keV jFinoguenov et al.l 
120011) . One nterpretation is again that metal loss due to starburst 
winds could be efficient in groups but that this process must have 
operated with near-uniform efficiency across the temperature range 
of our sample. While this may help to reduce ICM abundances in 
groups, it cannot, on its own, explain the SN II trends in Fig.jS] and 
much less so for the SN la products. On the basis of the available 
data, however, we cannot exclude the possibility that the discrep- 
ancy with respect to massive clusters could contain a contribution 
from this effect. 



4.4.2 AGN -driven winds 

The idea of central outflows driving material from group/cluster 
cores has received much recent attention with the observation that 
AGN in central cluster galaxies may be able to carve out cavi- 
ties or bubbles in the surro unding ICM through radio outbursts 
jMcNamara & NulsedfioOTl) . Evidence of such activity has also 
been identified in small groups, including several of the systems 
withi n our sample jAllen et al.j200^ : lMorita et al.l2006l . |jetha et al.l 
I2OO8I : iGastaldello et al.l 120091) . "ah indication that some AGN- 



induced gas mixing may have occurred in most, if not all, our 
groups is further provided by the radial extent of the central abun- 
dance excess, which is much broader than the distribution of optical 
light from the BGG (cf Rebuscg et al. 2006). Whether AGN-driven 
bubbles can rise far enough in clusters to displace l arge amounts of 
material from cluster cores remains unclear (e.g., iRoediger et al.l 
|2007|) . but semi-analytic models of galaxy formation have demon- 
strated that a variety of observational constraints in the group 
regime can be reasonably well satisfied in a scenario in which AGN 
may eject significant amounts of hot gas from low-mass haloes 
teower. McCarthv & Bensonir2008l) . If so, AGN would provide an 
attractive explanation for the gas and metal deficiency evident for 
such systems in Fig. [8] 

The standard assumption is that the observed cavities or bub- 
bles have arrived at their current location mainly by means of buoy- 
ancy forces. If bubbles of enriched gas have travelled from group 
cores to beyond rsoo by such forces only, the initial bubble entropy 
So = T/rig^^ must have been comparable to the ICM value Ssoo 
at rsoo • We can evaluate the temperature difference corresponding 
to AS = Ssoo — So and hence the energy required to heat the rel- 
evant gas mass in group cores for it to be buoyantly displaced be- 
yond rsoo. In a similar vein to the previous Section, we assume that 
the gas thus removed must account for the lower Jf-band IMLRia 
in the groups when compared to NGC6338. To treat all groups 
equally, we assume T500 = 0.67(r) for the ICM temperature at 
'"soo (cf. Paper I, equation 7), along with Z — Zq for the removed 
gas as observed in the group cores at present. For the initial entropy 
of the gas to be displaced, we note that inspection of our Chan- 
dra images of NGC741, HCG62, and NGC5044, all groups with 
known X-ray cavities, suggests an average projected distance from 
cavity centre to the optical centre of the BGG of roughly two-thirds 
of the D25 semi-major axis. Given our sample mean of ~ 22 kpc 
for the latter, this estimate is in excellent agreement with the mean 
projected distance of ~ 16 kpc seen for the radio-filled cavities 
in the c ombined g rou p and cluster sample of Blrzan et al. (200^1 
see also lRaffertv et aljbooeh . We therefore evaluate the initial gas 
density and entropy at r = 0.77325. 

Results are plotted in Fig. [18] in the form of the required 
heating energy AE. NGC4325 and NGC6338 itself are again ex- 
cluded, for the reason discussed earlier. Most groups scatter within 
a factor of a two in this plot, but NGC 2300 and 4125 are outliers 
due to their fairly high central entropies. The logarithmic mean of 
the heating energy corresponds to i5 ~ 8 x 10^" erg s^^, with no 
significant dependence on (T) {uk = —0.1) or Lk.bgg (o"a' = 
+ 1.5). While such energy requirements may excee d typical AGN 
outbu rst energies inferred for local clusters (e.g., iRaffertv et al.l 
l2006l) . the required energy could have been deposited in the course 
of multiple outbursts and so may be broadly consistent with ex- 
isting observational estimates. We note, however, that the associ- 
ated energy per total ICM particle is generally somewhat in ex- 
cess of typical values required by preheating scenarios (mean of 
AiJ ~ 2.4 keV particle"^), and does show some evidence for 
a negative trend with (T) {ok ~ —1.8) but not with Lk.bgg 
(CTK = 0.0). 

These results suggest that the energetic requirements for dis- 
placing the required amount of gas and metals by buoyancy are 
generally large for the current ICM configuration, but perhaps not 
prohibitively so. Another way to arrive at this conclusion comes 
from considering the implied efficiency e with which the central 
supermassive black hole (SMBH) must have converted accreted 
mass into outflow energy, assuming that the released energy rep- 
resents a converted rest mass corresponding to the current mass of 
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Figure 18. (a) Energy required to heat enriched gas in the group cores suf- 
ficiently for convection to carry this gas beyond rsoo and so achieve the 
same M-pa/ Lk ratio for SN la products within rsoo as in NGC6338. (b) 
The corresponding required black hole energy conversion efficiency for the 
estimated BGG black hole masses, (c) The corresponding time-averaged 
heating power {P) assuming an energy injection time of 10 Gyr (d) As (c), 
but compared to the estimated current mechanical power Lmoch of the cen- 
tral AGN, with typical constant (P) / Lmcch ratios marked by dashed lines. 
Uncertainties reflect the fractional errors on IMLRia. 



the SMBH. Estimating central SMBH masses from the BGG cen- 
tral stellar velocity dispersions (tq,* listed in Table|2l using 



M] 



1.2 X 10'' 



cro,. 



Mr. 



(5) 



200kms~^ 

jGebhardt et"ai]|200(]i) , we show the corresponding conversion effi- 
ciencies £bh = Ai?/MBHC^ in Fig. 118b. Results fall in the range 
e ~ 0. 1-5 per cent, with a median value of ~ 1 per cent. Interest- 
ingly, but perhaps coincidentally, this is in broad agreement with re- 
cent estimates of the actual SMBH efficiency in converti ng accreted 
rest mass into outflow energy , e ~ 0.5-2 per cent (e.g.. lAllen et al.l 
l200dlMerloni & Heinzll2008h . 

Further insight can be gained from comparing the required 
time-averaged heating power (P) to estimates of the current me- 
chanical AGN output in the groups. Fig.llSb shows the mean power 
associated with the derived AJ5 for a conservative assumption for 
the energy injection time-scale of 10 Gyr (again corresponding to 
an initial redshift z ~ 2). Fig. HSU compares this to rough esti- 
mates of the currently available power, based on the 1 .4 GHz radio 
power Pi. 4 of any central radio source as listed in Table[2] We con- 
verted Pi. 4, where available, into current mechanical AGN power 
Lmech ~ lO'^^L^alto, assumiug L.adio ~ (8.6 X 10'-* Hz)Pl.4 
as is appropriate f or a power-law radio spectrum of index unity 
( iBirzan et al1l2004) . Results reveal that current central AGN activ- 
ity in the groups is typically an order of magnitude weaker than the 
inferred time-averaged power required to accomplish the necessary 
gas removal. It is conceivable, however, that this may at least partly 
be related to selection effects, since any system with a powerful ac- 
tive radio source strongly disturbing the ICM might not have made 
it into our sample, which only includes groups which appear rea- 
sonably undisturbed on large scales in Chandra data. 



While acknowledging that the estimates in Fig.[T8]result from 
rather simplistic considerations which, among other things, ignore 
the potential contribution from AGN not associated with the BGG, 
the derived energy requirements suggest that some buoyant dis- 
placement of gas to beyond rsoo by AGN activity may have taken 
place in the recent past, but the majority of any such activity is 
likely to have occurred when the ICM entropy configuration was 
more conducive to blow-out. Also note that our energy estimates 
are based on total mass losses assuming solar abundances for the re- 
moved gas, but that larger amounts of gas (and hence lower metal- 
licities and larger outflow energies) would be needed if AGN were 
to also explain the general shortfall of hot gas compared to clusters 
evidenced by Fig. [13] In combination, this argues in favour of gas 
and metal loss predominantly, but not necessarily exclusively, at 
earlier epochs, closer to the redshift peak of metal production and 
possibly prior to group collapse. 

We do note, however, that recent AGN outflows could eas- 
ily have displaced highly enriched gas from group cores and so 
have contributed to the central Z-T trend in Fig. [TO] Using sim- 
ilar assumptions as those underlying Fig. [TS] but only considering 
the 'missing' Fe mass within O.lrsoo as derived from the data in 
Fig. 1 12b. and assuming that this gas has simply been expelled be- 
yond O.lrsoo, the resulting energy requirements display sample 
mean values of (logA_E) « 59.6 erg, (r^BH) ~ 0.1 per cent, 
(logP> ^ 42.1 erg and ((logP)/(logPn,ech)> ~ 1.0. These 
requirements appear perfectly reasonable, and suggest that current 
AGN power in the groups is generally sufficient to buoyantly dis- 
place the required amount of enriched gas from within O.lrsoo. 
The requirements are generally slightly harder to satisfy in hotter 
systems (as required by FigllOll, with both required logAiJ and 
(log P)/(log Pjnoch) showing a weakly significant tendency to in- 
crease with T, at crx = +1.4 and cjk = +1.6, respectively. 



5 PREHEATING AND METAL EJECTION FROM 
FEEDER FILAMENTS 

The discussion in the preceding Section focussed on explaining 
the low IMLR in groups in the framework of a decreased poten- 
tial for in situ metal production/release or retention in present-day 
group-sized haloes. Motivated by the limited success of these ex- 
planations, with the possible exception of AGN feedback, we now 
consider the possibility that it is mainly related to processes occur- 
ring at the epoch of galaxy formation. For example, the AGN re- 
sults in Section |4!4l demonstrate that substantial, although not nec- 
essarily prohibitive, amounts of energy need to be injected into the 
ICM to accomplish the necessary gas and metal removal at low 
redshift from these groups. It would, however, have been energeti- 
cally favourable to remove some of this material at a much earlier 
stage, when potential wells were shallower, and the gas more ten- 
uous and likely exhibiting a more irregular geometry much more 
conducive to gas removal. It seems attractive to connect this pos- 
sibility with the result that both star formation and AGN activity 
were much stronger at high redshift, both peaking around z ~ 2-3 
(Bouwens et al. 2007; Hopkins, Richards & Hernquist 2007). This 
section summarizes evidence that a significant fraction of metals in 
groups must have been generated before or during this phase, and 
discusses to what extent proto-galactic winds may have heated and 
dispersed this material beyond what is now rsoo. 
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5.1 The case for early enrichment 

There is substantial evidence that a significant fraction of the 
present-day metals in the Universe must have been generated at 
very early epochs. Enrichment of the intergalactic medium was tak- 
ing place already at 2 ~ 6 tSimcoe 2006 ). with a non-negligible 
component of the produced metals residing in Lyman break galax- 
ies, some of which presumably evolved into the massive cluster 
ellipticals seen today. Contenders for this early enrichment include 
proto-galactic winds occurring when (what are now) massive early- 
type galaxies formed the bulk of their stellar mass. In galaxy clus- 
ters, evidence for substantial enrichment prior to cluster collapse 
has b een reported on the basis of measured large-radius abundances 
jpuii ta et al. 2008). More generally, results for large cluster sam- 
ples indicate that the hot gas had achieved a mean Fe abundance of 
Zpe « 0.3 Z0 by 2 ~ 1, comparable to that seen in nearby clus- 
ters, althou gh some evolut ion in the central Fe content has occurred 
since then ( Balestra et al. I l2007l : iMaughan et al."2008V The latter 
result draws support from the simulations of Cora et al. ( 2008), 
which suggest that cluster ICM abundances at large radii (~ rsoo) 
were largely established at z > 1, with only central abundances in- 
creasing significantl y since then. In small er systems, the situation 
may be less clear (Finog uenov et alj|2002h but in the specific case 
of our groups, a couple of results do suggest that the material in 
their outskirts was predominantly enriched at early epochs by core- 
collapse SNe, whereas the abundance excesses in group cores may 
have been generated more recently. 

First, the inferred number ratio of SN la to core-collapse SN 
in group cores of ^ 0.4 is consistent with that measured in the lo- 
cal {z < 0.7) Universe in the Hubble and Chandra Deep Fields, 
but disagrees with predictions at z > 1 for reasonable evolutionary 
models (Dahlen et a l. 2004). In the absence of serious recent gas 
mixing (supported by the relatively undisturbed X-ray morphology 
and the pronounced Fe gradients of our groups), this suggests that 
the SN la-dominated group cores in general, and the central Fe and 
Si excesses in particular, have been generated within the last ~5- 
7 Gyr The abundance pattern seen in the heavily SN Il-dominated 
outskirts which represent most of the ICM is, on the other hand, 
consiste nt with expe ctations for the SN la/SN II ratio at redshifts 
z > 1 toahlen et a i. 2004). Hence, most of the SN II activity in 
these groups likely took place at high redshift. This is in accordance 
with the result that most of our groups show a galaxy morphology- 
dens ity relation of comparable strength to that observed in clus- 
ters faelsdon & Ponmanl l2003h. suggesting that star formation and 
hence SN n activity has largely ceased at present in many of the 
bright galaxies in the central regions. 

Second, even when allowing for enrichment through stellar 
mass loss from the central early-type in each group, the discussion 
in Section 13.31 shows that an abundance contribution from SN II 
is, on average, required everywhere inside rsoo. The near-uniform 
distribution of SN II products revealed by Fig. [6] suggests either an 
earlier enrichment epoch for these than for the SN la ejecta, allow- 
ing time for substantial mixing before and during group collapse, 
or a much more efficien t redistribution mechanism for SN II ejecta. 
However, the results of iFerrara, Pettini & Shchekinovl ([2000) sug- 
gest that the distribution of SN II products in the absence of 
other mixing mechanisms should roughly follow that of the galaxy 
light. While this is observe d for the central SN la products in 
at least some of our groups jMendel et al]|2009l) . the SN II met- 
als seem to trace more closely the much more extended distribu- 
tion of the hot gas. The idea that SN II enrichment largely oc- 
curred before group/cluster collapse, as advocated by, for exam- 



ple, |Finoguenov|et^] i200lh . implies that the metals seen in group 
outskirts must have been generated by SN II before a significant 
number of SN la began exploding. Depending on the exact delay 
time distribution of SN la, this would imply SN II iron enrichment 
to a level of ~ 0.1 Z© (Fig.|6l at redshifts z > 2 — 3. As discussed 
earlier, the gentle central rise seen for SN II products can then be 
explained if a relatively small fraction of these metals were released 
after the group collapsed, possibly facilitated by galaxy-galaxy and 
galaxy-ICM interactions, and with an additional contribution asso- 
ciated with stellar wind loss from the central early-type galaxy. 

By considering the typical time-scales for metal production in 
the groups, an additional argument can be provided in support of a 
prominent high-redshift contribution to ICM enrichment. For this, 
we focus first on the Fe produced by SN la, which is less likely 
than SN II products to become locked up in stars in the early-type 
galaxies that dominate the .ft'-band output in our groups. Using the 
same general approach as for Fig. [T2I while assuming a constant 
SN la rate per unit Lk equivalen t to the value observed in local 
early-types jMannucci et aljfioOS) . we plot the resulting enrich- 
ment times corresponding to the total if -band IMLRia within rsoo 
in Fig. 1 19b. Under these assumptions, the enrichment times should 
again be considered lower limits, because no account is made of 
metals locked in stars, metals that have left the group potential, 
or the growth of stellar mass within the group due to galaxy as- 
sembly and the addition of galaxies to the group over the time- 
scales involved. With results suggesting time-scales of > 10 Gyr in 
40 per cent of our groups, typically in the hotter and more SN la- 
dominated ones, it is far from clear that specific SN la rates seen in 
present-day early-types can account for the observed metal masses. 
The discussion in previous Sections further indicates that contribu- 
tions from stellar mass loss and intracluster stars are insufficient for 
bringing these time-scales below ~ 10 Gyr in all groups; either the 
specific SN la rate must increase with redshift, or 'external' enrich- 
ment must be invoked to explain the observed amount of Fe from 
SN la. 

This issue is even more acute for SN II, whose contribution 
to Fe enrichment is globally at least as important as that of SN la 
(cf. Table[T](, but for which onl y upper limits to the ob served rate in 
local early-types are available jMannucci et al.ll2005[) . These gen- 
erally imply enrichment times far in excess of the age of the Uni- 
verse. However, some morphological transformations and trunca- 
tion of star formation has l ikely taken place in the group members 
( iHelsdon & Ponmanll2003h . so assuming current early-type SN II 
rates could well underestimate the metal output over cosmic time. 
In a crude attempt to correct for this, we employ the conservative 
assumption that all stellar mass not currently in the BGG has pro- 
duced SN II according to the rate in local late-type spirals. The 
result is shown in Fig. 1 19b . As is clear, time-scales are still pro- 
hibitively large in several cases, implying that SN II rates must also 
evolve with redshift, and even more so than for SN la. Specific 
star formation rates, and hence metal production rates, must have 
been considerably higher in the past in these groups. Though not 
a surprising result, it lends independent support to the notion that 
an important fraction of ICM metals within rsoo was generated at 
early epochs. 

5.2 Linking enrichment to preheating 

With much of the ICM pollution likely taking place around z ~ 2- 
3, close to the redshift peaks of the cosmic star formation rate and 
of the AGN luminosity density, it is conceivable that metals and 
energy were released from galaxies in substantial amounts before 
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Figure 19. (a) SN la enrichment times based on the X-band IMLRia. plot- 
ted against the fractional Fe mass from SN la. (b) Similar for SN II, as- 
suming local Sbc/d SN II rates for all satellite galaxies, as a function of the 
BGG contribution to total K-hand light. 



gas and galaxies collapsed into groups, i.e. while still located in fil- 
aments of low overdensities. This would not only suggest that the 
trends in Fig. [8] were at least partly established at this epoch, but 
also offer a possible explanation of their origin in gas and metal 
loss from the filaments that fed the growth of the group structures 
observed today. Accretion along filaments implies, unlike the case 
for virialized structures, a preferred direction along which the ejec- 
tion and loss of galactic material is strongly favoured. Given the 
available data, we can only make some general speculations about 
this possibility, but it seems an attractive one in light of the some- 
what limited success of the other explanations considered. 

ICM enrichment at high redshift could possibly be linked to 
the putative 'preheating' invoked in some models, by which the 
ICM in these systems is assumed to have had its entropy raised 
at high redshift in order to explain the similarity breaking in the 
X-ray scaling properties of groups and clusters. Simulations and 
observations (Borgani et al. 2008a and references therein) indicate 
that the associated energy required is of order ~ 1 keV per parti- 
cle, depending on the exact redshift of energy injection. It cannot 
greatly exceed this value, as preheating would otherwise have de- 
stroyed the ICM in T < 1 keV groups. While the initial argument 
for introducing this was to explain the seemingly elevated entropy 
levels in group cores, preheating could also affect the gas at larger 
radii. For example, the low gas mass-to-light ratio within rsoo in 
cool systems could be related to the possibility that groups have se- 
riously inflated ICM distributions as a consequence of preheating 



through early galactic winds. The resulting entropy of the heated 
(and possibly enriched) gas would make it more resilient to gravi- 
tational compression, to a certain extent preventing its re-collapse 
into the group structures seen at low redshift. Thi s is somewhat akin 
to the 'preferential infall' scenario proposed bv iFinoguenov et al] 
in which small groups preferentially accrete low-entropy, 
low-metallicity gas that was only mildly heated at high redshift. 

One implication would be that a relatively larger fraction of 
any gravitationally bound X-ray gas and metals should reside be- 
yond rsoo in cooler systems. We are not in an optimal position to 
test this, as it would require sensitive X-ray data extending beyond 
this radius for all groups. We can only note that with the adopted 
density parametrizations, the ratio of gas mass within O.Srgoo to 
that within 0.5 — Irsoo actually shows a weak anft-correlation 
with (T) ((J A' = —1.6) for our groups, consistent with the cor- 
relation seen between /3 and T for the GEMS sample in general 
(o — —1.4; Osmond & Ponman 2004), from which the majority of 
our groups were selected. However, iSanderson & PonmanI ( l2003h 
do find a positive dependence of /3 on T when including clusters in 
the comparison. The group-cluster discrepancy within 7-500 could 
contain some contribution from this effect, and if so would sug- 
gest that a more fair comparison of systems across a wide range of 
masses should be conducted within radii of lower overdensity such 
as r2oo- 

Our speculation that a considerable fraction of the (centrally 
concentrated) SN la products were released after group collapse 
suggests that at least part of the enrichment is not associated with 
preheating but must have occ urred at later time s at z < 1. This 
agrees with the conclusion of IFinoguenov et al] ( l2002h . based on 
the absence of a correlation between entropy at radii of equal en- 
closed mass vs. Fe or Si abundance across a wide range of system 
masses. However, SN II provide the clear majority of the estimated 
SN energy released in our groups, with an average contribution to 
EsN,* in Fig. [15] of ~ 0.8 keV per particle, most of which was 
presumably released in early starbursts prior to group collapse. A 
potential caveat to the assumption that this value also corresponds 
to the energy delivered to the ICM by galaxies is that much of the 
SN explosion energy may have been rapidly radiated away in the 
initial starburst phase rather than released from galaxies in kinetic 
form. Whereas radiative losses are presumed to be minor at low 
redshift, with direct escape of the SN explosion energy facilitated 
by the presence of bubbles and chimneys in the ISM carved out by 
previous generations of SN (and winds from their progenitors), the 
ISM was denser, colder, and less porous at high redshift, possibly 
enhancing radiative losses. If so, the ~ 0.8 keV per particle inferred 
from the stellar mass of the group members may represent an upper 
limit to the actual SN II energy deliv ered to the galaxy surround- 
ings. On the other hand, simulations jStrickland & Stevensllioooh 
indicate that coordinated SN explosions, such as those likely tak- 
ing place during vigorous high-redshift starbursts, tend to lower 
radiative losses. Furthermore, Fig. [15] suggests that the estimated 
SN energy associated with the ICM metal mass broadly matches 
the 'available' energy, so any overestimation of the SN energy im- 
parted to the ICM should be modest. 

Even if only a minor fraction of this energy were actually de- 
livered to the surroundings in proto-galactic winds, it could still 
have contributed substantially to pre-heating the filaments feeding 
group formation. It may not be sufficient to entirely explain the sim- 
ilarity breaking between groups and clusters, or the presence of an 
'entropy ramp' in group and cluster cores, but it certainly eases the 
energetic requirements imposed on any other source to deliver the 
remaining heating, and it may leave limited scope for AGN feed- 
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back. In the cluster regime, numerical simulations jTomatore et al.l 
|2004) suggest that a feedback mechanism not associated with SN 
has to be acting at early stages of cluster formation, before most of 
the mass is accreted in the cluster potential well, in order to avoid 
the centrally peaked Fe distribution seen in some simulations. At 
lower masses, however, where our results show stronger Fe gra- 
dients that in clusters, the need f or AGN feedback is perhaps less 
pressing. The simulation work of lOave et al. 1 1I2OO8'), for example, 
indicates that elevated entropy levels in groups may be largely ex- 
plained by radiative cooling combined with SN-driven outflows, 
with no strong need for AGN-driven preheating. 

Regardless of the heating source, the energetic requirements 
for blow-out of gas via buoyancy discussed in Section l4!4l would 
be strongly relaxed if the surrounding material is at a lower over- 
density than currently seen in group cores. For example, Fig.llSlis 
based on assuming current ICM densities within ~ 10 — 20 kpc of 
the group cores, which generally translate into cosmic overdensities 
of order 10^. If the gas to be heated has lower overdensities by 
a factor ~ 100 (certainly realistic for pre-virialized structures such 
as inflowing filaments), the requirements for heating it to current 
entropies at rsoo would be reduced by a factor of ~ 20, i.e. to just 
~ 0.1 keV particle"^ on average. This should be within range of 
either SN or AGN at 2: ~ 2 — 3, suggesting that the energy required 
to drive metals out of filaments is, indeed, available. 

Whether the balance between bulk gas loss and selective 
loss of both SN la and SN II products suggested by our results 
can be achieved is a non-trivial question though. The shortfall of 
SN II products in groups seems largely explicable on the basis 
of pre-collapse metal loss in starburst winds, but removal of large 
amounts of Fe from SN la at early epochs may also require some 
of this material to have been synthesized very rapidly, for exam- 
ple by a population of SN la exploding promptly after star for- 
mation. The existence of such a population has been proposed by 
iMannucci, Delia Valle & Panagia (2006) on the basis of observed 
SN la rates at high redshift, but the ne ed for this popula tion at 
early times is still debated (Dahlen, Strol ger & Riessll2008t) . Nev- 
ertheless, even in the 'standard' SN la delay-time scenario, a non- 
negligible fraction of SN la products may have been released at 
fairly high redshift. For example, the model of .Cor a ( 2006) pre- 
dicts a peak in the SN la iron injection rate at z ~ 2 — 3 in satellite 
galaxies and even earlier in the central galaxy. 

While some SN la material may thus have been ejected from 
cool systems in early starburst winds, our discussion of Fig. [6] 
which shows that the SN la contribution becomes insignificant out- 
side ~ O.Trsoo still argues against a significant pre-collapse SN la 
contribution and so against substantial loss of SN la metals from fil- 
aments. To explore this issue further, we plot in Fig.|20]Fe masses 
from SN la within different radial intervals, all simply normalised 
by total stellar light Lk ■ The plot suggests that, over the mass range 
considered here, the reduction in IMLRia in lower-mass systems 
has preferentially occurred within the central group regions, with 
no trend with (T) seen beyond O.Srsoo- This preferential loss of 
SN la products from the central Fe peak seems consistent with post- 
collapse AGN outflows (and possibly some cooling) as important 
contributors to the shortfall of SN la products in cool systems. The 
considerations in Section [4.4. 21 then suggest that this loss cannot 
generally have occurred very recently but rather must have taken 
place over considerable time-scales following group collapse. 
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Figure 20. Mass of Fe from SN la within different radial ranges, normalised 
by total Lif. As in Fig. 1 101 only groups with data covering the full radial 
interval are shown in each panel. Labels in lower right comers specify the 
significance aj^ of any linear correlation. 



6 DISCUSSION AND CONCLUSIONS 

Based on the results pr esented in Paper I in this series 
dRasmussen & Ponmanl2007>) . in this second of two papers we have 
investigated the distribution and total mass of SN la and SN II prod- 
ucts in the ICM of 15 groups, for comparison to the total ICM mass 
and to the optical luminosities and stellar masses of the group mem- 
bers. We re-emphasize here that our group sample was not intended 
to be statistically representative in any sense. Nevertheless, the re- 
sults should still provide useful constraints on models of the evo- 
lution of galactic systems on scales from small groups to massive 
clusters. We note that there are features of the observed abundance 
distributions that we have not covered in any detail. These include 
the broad extent of the central abundance excess, reaching well be- 
yond the distribution of stellar light from the central galaxy, as well 
as the curious abundance dips seen for both Fe and Si at the very 
centre for several groups, manifest also in the averaged profiles in 
Fig. [T] Any detailed discussion of these features is likely to bene- 
fit from comparison of group properties to those of non-cool core 
systems, however, and so is deferred to a future paper. Our existing 
results can be summarized as follows: 

1. As anticipated, the total mass of iron and silicon in the hot 
gas within rsoo scales strongly with total group mass as measured 
by the mean X-ray temperature outside any cool core. Values range 
fromMpc = (0.4-23) X 10* MoandMsi = (0.6-35) x 10* Mq 
across the 0.3-2.1 keV temperature range of the sample. SN II have 
provided on average ~ 60 per cent of the Fe and ~ 90 per cent of 
the Si inside rsoo, with a mildly significant tendency (at ^ 2(t) 
towards increasing SN II dominance in cooler systems. 

2. For the sample as a whole, the distribution of SN la products 
is strongly peaked towards the group core, whereas SN II enrich- 
ment is nearly uniform throughout the ICM. The Fe mass within 
O.lrsoo (roughly the size of the cool core, where present) can be 
explained by prolonged enrichment by SN la in the central early- 
type galaxy, aided by stellar mass loss and intracluster stars. In par- 
ticular, if adopting a plausible model for the cosmic evolution of 
the SN la rate, typical enrichment time-scales are < 5-7 Gyr in all 
groups, comparable to similar results for central cluster galaxies. 
The mild central rise observed for SN II products can be explained 
by stellar mass loss from the central galaxy, potentially with an ad- 
ditional contribution from metals released by other group galaxies 
via density-dependent release mechanisms such as galaxy-ICM in- 
teractions. 

3. Due to the combination of lower gas mass-to-light ratios 
and metal abundances in cooler systems, resulting aggregate metal 
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mass-to-light ratios within rsoo show an order-of-magnitude varia- 
tion across the sample, only approaching results for massive clus- 
ters at the high-temperature end. For a given amount of stellar light, 
there is thus a substantial shortfall of X-ray metals in cooler sys- 
tems, a deficiency which is found to be present both for SN la 
and SN II products. This suggests that metals from both SN types 
were either synthesized or released from galaxies less efficiently 
in lower-mass systems, that metals were lost from such systems, or 
that the intergalactic gas accreted by small groups was less enriched 
than that accreted by more massive systems 

4. Both Fe and Si abundances in group cores (r < O.Srsoo) 
show a positive trend with group temperature which is significant at 
~ 3(7. No corresponding trend is seen at larger radii (r > O.Srsoo) 
however, arguing against recent preferential loss of enriched mate- 
rial from cooler systems as an explanation of their lower global 
metal-mass-to-light ratios. Estimates of the current mechanical 
power of the central AGN in the groups indicate that the metal- 
licity trend in group cores may at least be partly due to the expul- 
sion of enriched material to beyond O.lrsoo by AGN-driven buoy- 
ant bubbles, perhaps with some contribution from radiative cooling 
of enriched gas in the core, causing less enriched material to flow 
inwards. Alternatively, it can be understood as due to the action 
of a metal release mechanism whose efficiency depends on local 
ICM density and the total mass of the system, such as ram pressure 
stripping, but it cannot be explained by a systematic variation in 
the contribution from intracluster stars with group temperature. The 
absence of a similar reported trend in systems above T ~ 2 keV 
may indicate that the relevant processes 'saturate' in cluster cores, 
perhaps as a consequence of near-complete gas removal from core 
galaxies in such environments, or due to the larger heating energies 
or time-scales required to accomplish significant metal removal by 
AGN outflows or radiative cooling in such systems. 

5. As a consequence of the lower Fe abundances and mildly 
increased contribution from SN II in cooler systems (i.e. higher 
a/Fe ratios), we find tentative evidence for the ICM within the 
central few hundred kpc (r < O.Srsoo) to appear chemically 
less 'evolved' in these groups. This may reflect a more extended 
star formation history in less massive groups, consistent with the 
slightly bluer average galaxy colour seen in these systems. 

6. Metal release from galaxies to the ICM appears to have 
been efficient in these groups, with, on average, an estimated ~ 
70 per cent of the metals synthesized in galaxies currently present 
in the ICM. While thi s is in fair agreement wit h recent independent 
estimates in clusters feivanan dam et alj |2009^. there is also some 
evidence for the release fraction to be systematically lower, and 
well below cluster values, in the cooler groups in our sample. 

7. Among possible explanations for the metal deficiency in 
groups within rsoo, we can rule out radiative cooling as an impor- 
tant gas and metal sink. Some contribution may instead come from 
the more extended star formation histories and lower metal release 
efficiencies in cooler systems, as well as the tendency for such sys- 
tems to contain a relatively larger fraction of their gas beyond rsoo • 
Within the current ICM configuration, the energetic requirements 
for the necessary gas and metal loss induced by AGN- or SN-driven 
winds are rather large (although not necessarily prohibitively so), 
generally favouring loss of material at the incipient stages of group 
evolution. Other mechanisms such as thermal conduction transport- 
ing gas and metals out beyond rsoo in cooler systems can probably 
be ruled out, due to the strong temperature dependence anticipated 
for heat conduction unsuppressed by magnetic fields (Spitzer con- 
ductivity oc T^^^). 

8. Comparison to observed ratios of SN la to core-collapse SN 



in deep field data jPahlen et all2004l) shows that group cores dis- 
play abundance ratios consistent with inferred SN ratios in the low- 
redshift Universe {z < 0.7) but are inconsistent with predictions 
for 2 > 1. This argues in favour of relatively recent enrichment 
of the central group regions, including the build-up of the central 
abundance excesses over the past < 6-8 Gyr In contrast, the heav- 
ily SN Il-dominated abundance ratios in group outskirts does agree 
with predictions for z > 1, suggesting that most of the ICM was 
enriched to ~ 0.1 solar at redshifts close to the peak of the cosmic 
star formation rate. The near-uniform radial distribution of SN II 
products corroborates this picture by arguing in favour of substan- 
tial SN II enrichment prior to or during group collapse. 

9. This raises the possibility that powerful starburst winds at 
high redshift have contributed significantly to the gas and metal 
deficiency in groups, causing loss of baryons from the filaments 
feeding the growth of these structures. The estimated SN explosion 
energy associated with high-z star formation in the group mem- 
bers seems sufficient for this, and may further have provided a sub- 
stantial fraction of the preheating energy required by some models 
for explaining the X-ray scaling properties of groups. If so, it is 
not clear that a prominent contribution to preheating from quasar 
activity is required in these systems. Post-collapse AGN feedback 
can nevertheless still have played a prominent role in subsequently 
ejecting mainly SN la metals from group potentials. 

One implication of our results is that present-day clusters can- 
not have formed through hierarchical build-up of smaller struc- 
tures resembling the groups studied here. This is, of course, already 
clear from the well-established differences in cold fraction (cf. Sec- 
tion |4. It , but the deficiency of gas-phase metals in smaller systems 
poses an additional problem for this scenario. Of course there is 
no compelling reason to believe that present-day groups are chem- 
ically similar to those that participated in cluster build-up at higher 
redshift, although our considerations suggest that the metal defi- 
ciency in low-mass systems has not been established recently. Al- 
ternatively, the proto-clusters into which such groups merged have 
been massive enough to retain or re-accrete enriched material shed 
by their smaller progenitors. Deep X-ray observations of high-2; 
groups may eventually shed light on these issues, but the questions 
are unlikely to be settled by the current generation of X-ray satel- 
lites. Numerical work may also provide clues to the chemical prop- 
erties of cluster progenitors and in which ways these depart from 
those of our groups. A prerequisite for such a comparison, how- 
ever, is that models can accurately reproduce observational results 
for the ICM abundance distribution in groups at low redshift. While 
significant improvements have been made in the past few years in 
this regard, some issues remain such as reprodu cing the observed 
distribution of a-elements (see, e.g.. lDave et al.l l2008). 

Alternatively, our groups are conceivably not representative 
of (X-ray) groups in general. Although exhibiting X-ray luminosi - 
ties that span two orders of magnitude jOsmond & Ponmarill2004l) . 
our group sample is still largely restricted to cool-core groups that 
appear morphologically undisturbed. In order to gain a more com- 
plete picture of the enrichment and feedback history in groups, sim- 
ilar analyses would need to be extended to a more heterogeneous 
group sample, including a substantial number of systems without a 
cool core. At higher masses, such systems tend to show much flatter 
abundance gradients than those reported here (e.g., De Grandi et al] 
2004), suggesting fundamental differences in ICM enrichment and 
heating histories between cool-core and non-cool core systems. 
Such a comparison is the subject of ongoing work (Johnson et al., 
in preparation), and may ultimately help in providing an improved 
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understanding of the chemical and thermodynamical evolution of 
baryons in both galaxies and the intergalactic medium. 
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